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ABSTRACT 

We report results from a large molecular line survey of Luminous Infrared Galaxies (LIRGs: 
\ LiR>10 n L0) in the local Universe (z^O.l), conducted during the last decade with the James 

Clerk Maxwell Telescope (JCMT) and the IRAM 30-m telescope. This work presents the CO 
and 13 CO line data for 36 galaxies, further augmented by multi-J total CO line luminosities 
available for other IR-bright galaxies from the literature. This yields a combined sample of 
psj ■ N=70 galaxies with the star-formation (SF) powered fraction of their IR luminosities span- 

ning l|^'^(10 10 -2x 10 12 ) L-0, and a wide range of morphologies. Simple comparisons of 
their available CO Spectral Line Energy Distributions (SLEDs) with local ones, as well as 
radiative transfer models discern a surprisingly wide range of average ISM conditions, with 
most of the surprises found in the high-excitation regime. These take the form of global CO 
SLEDs dominated by a very warm (Tkm^lOOK) and dense (n^lO 4 cm -3 ) gas phase, involv- 
O'N ing galaxy-sized (^(few)x 10 9 M Q ) gas mass reservoirs under conditions that are typically 

found only for ^(1-3)% of mass per typical SF molecular cloud in the Galaxy. Furthermore 
some of the highest excitation CO SLEDs are found in Ultra Luminous Infrared Galaxies 
(ULIRGs, Lir>10 12 Lq), and surpass even those found solely in compact (star-formation)- 
powered hot-spots in Galactic molecular clouds. Strong supersonic turbulence and high cos- 
mic ray (CR) energy densities rather than far-UV/optical photons or SNR-induced shocks 
r> . from individual SF sites can globally warm the large amounts of dense gas found in these 

merger-driven starbursts and easily power their extraordinary CO line excitation. This excit- 
ing possibility can now be systematically investigated with Herschel and ALMA. As expected 
for an IR-selected (and thus (SF rate)-selected) galaxy sample, only few "cold" CO SLEDs 
are found, and for fewer still a cold low/moderate-density and gravitationally bound state 
(i.e. Galactic-type) emerges as the most likely one. The rest remain compatible with a warm 
and gravitationally unbound low-density phase often found in ULIRGs. Such degeneracies, 
prominent when only the low-J SLED segment (J=l-0, 2-1, 3-2) is available, advise against 
using its CO line ratios and the so-called X CO =M(H2)/L CO (1-0) factor as star-formation mode 
indicators, a practice that may have lead to misclasification of the ISM environments of IR- 
selected gas-rich disks in the distant Universe. Finally we expect that the wide range of ISM 
conditions found among LIRGs will strongly impact the X co factor, an issue we examine in 
detail in paper II (Papadopoulos et al. 2012). 

Key words: galaxies: ISM - galaxies: starburst - galaxies: AGN - galaxies: infrared - ISM: 
molecules - ISM: CO 



1 INTRODUCTION 
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The population of luminous infrared galaxies (LIRGs), discovered 
by the Infrared Astronomical Satellite (IRAS) to have bolometric 
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luminosities dominated by the infrared part of their Spectral En- 
ergy Distributions (SEDs) (e.g. Soifer et al. 1987), contains some 
of the most extreme star-forming systems in the local Universe. At 
Lir > 10 11 L,0 these deeply dust enshrouded star-forming sys- 
tems dominate the luminosity function of the local Universe, and 
at Lir ^ 10 12 Lq surpass even optically selected QSOs (Soifer & 
Neugebauer 1991; Sanders & Ishida 2004 and references therein). 
Their large reservoirs of molecular gas mass (~1O 9 -1O 1O M ) dis- 
covered via CO J=l-0 observations (Tinney et al. 1990; Sanders, 
Scoville, & Soifer 1991; Solomon et al. 1997), along with clear evi- 
dence of strong dynamical interactions and mergers in many LIRGs 
(e.g. Sanders & Ishida 2004), make these systems unique local ex- 
amples of dust-enshrouded galaxy formation in the distant Universe 
(e.g. Smail, Ivison, & Blain 1997; Hughes et al. 1998). 

Interferometric imaging of the CO 1-0 line (and occasionally 
of J=2-l), revealed gas disks with D<0.5kpc and surface densi- 
ties £(H 2 )~(10 2 -10 4 ) M Q pc~ 2 , often decoupled from the stellar 
components of merging galaxies (e.g. Sanders et al. 1988a; Wang 
et al. 1991; Planesas et al. 1991; Bryant & Scoville 1996, 1999; 
Downes & Solomon 1998). Naturally the CO(2-1)/(1-0) ratio was 
the first to be systematically measured (Krugel et al. 1990; Braine 
& Combes 1992; Horellou et al 1995; Aalto et a. 1995; Albrecht, 
Krugel, & Chini 2007) even if it is insensitive to the presence of 
dense and warm gas typical near SF sites. The advent of submm 
interferometric CO(3-2) imaging revealed the distributions of such 
gas in LIRGs via the distribution of CO (3-2)/(l-0), (3-2)/(2-l) 
ratios (e.g. Sakamoto et al. 2008; Wilson et al. 2009; Iono et al. 
2007, 2009), but its utility remains limited by dissimilar u-v cover- 
age and lack of zero-spacing information (Iono et al. 2004). Thus 
single dish measurements of total molecular line luminosities and 
their ratios remain a primary tool for probing the average molecu- 
lar gas properties in LIRGs (and a prerequisite for interferometric 
images that contain all spatial information). 

Ideally a combination of low to mid-J rotational lines of 
heavy rotor molecules with high critical densities such as HCN 
(n cr ~(2-40)xl0 5 cm~ 3 for J=l-0, 3-2), and CO J+HJ lines 
from J=0, 1, up to at least J^2 (E 3 /k B ~33 K, n crit ~10 4 cm" 3 ) 
are necessary to probe the large range of physical properties within 
GMCs (T kin ~(15-100)K, n(H 2 )~(few)x(10 2 -10 6 )cm _3 ). Sen- 
sitivity limitations and/or lack of multi-beam receivers confined 
such measurements to few nearby LIRGs and mostly towards their 
nuclei (e.g. Devereux et al. 1994; Dumke et al. 2001; Meier et al. 
2001; Zhu et al. 2003), while heavy rotor molecular lines are faint 
with HCN( l-0)~( 1/5- 1/10) x CO( 1-0) even in ULIRGs, and as low 
as ~ l/40xCO(l-0) in typical spirals (Solomon et al. 1992). Multi- 
J observations of the luminous CO line emission are not limited by 
sensitivity rather by the lack of multi-beam receivers and/or beam- 
matched observations in widely different frequencies. Thus, while 
several studies use CO(3-2)/(l-0) as a warm/dense gas tracer in 
substantial LIRG samples (e.g. Mauersberger et al. 1999; Yao et al. 
2003; Narayanan 2005; Mao et al. 2011), these are typically con- 
fined solely within their nuclear regions. CO J=4-3 or higher-J line 
observations are even more sporadic, hindered by increasing atmo- 
spheric absorptions at u> 460 GHz, while also confined to the nu- 
clear regions of nearby LIRGs (e.g. White et al. 1994; Glisten et al. 
1996; Petitpas & Wilson 1998; Nieten et al. 1999; Mao et al. 2000; 
Bayet et al. 2006). This is unfortunate since, while the difficulty 
of measuring reliable CO line ratios increases with wider J-level 
separations (for the same dish: f2j^™j/f2^Q m '=(J + 1)~ 2 ), so 
does their diagnostic power. Nevertheless molecular Spectral Line 
Energy Distributions (SLEDs) remain the key tool for probing the 
average state of the molecular gas in galaxies, and for estimating to- 



tal and star-forming molecular gas masses. Finally local CO SLEDs 
provide a necessary benchmark for the usually more sparsely sam- 
pled ones for galaxies at high redshifts (e.g. Weiss et al 2007), and 
have been recently used as SF mode (merger versus disk-driven 
star formation) indicators for gas-rich near-IR selected disks in the 
distant Universe (Dannerbauer et al. 2009; Daddi et al. 2010). 

We used the James Clerk Maxwell Telescope (JCMlEl) on 
Mauna Kea in Hawaii (USA), and the IRAM 30-meter telescope 
at Pico Veleta (Spain) to conduct a multi-J CO, HCN line survey 
of such systems and study the molecular ISM in some of the most 
extreme star-forming systems in the local Universe. First results re- 
garding the CO J=6-5 lines and dust emission SEDs have already 
been published (Papadopoulos et al. 2010a,b, hereafter P10a,b). In 
this paper we present the entire (JCMT)+(IRAM) CO, 13 CO line 
dataset, further augmented with all reliable measurements of to- 
tal CO line luminosities of LIRGs from the literature, which yields 
currently the largest such database assembled for local star-forming 
galaxies. For a few LIRGs in our sample their CO SLEDs have been 
extended towards higher-J levels (from J=5^1 up to J=13— 12) with 
the SPIRE/FTS aboard the Herschel Space Observatory. 

The layout of this work is as follows: in Sections 2, 3 we de- 
scribe the sample and the observations, Section 4 contains the data 
reduction, the literature search, and line flux rectifications for any 
discrepant values found. In Section 5 we investigate the molecu- 
lar ISM exitation in LIRGs using comparisons with local environ- 
ments, and propose a new two-phase ISM model for extreme star- 
bursts, motivated by recent views on SF-powered radiative feed- 
back onto the ISM. In the same saction we also briefly discuss the 
role of AGN in driving the global ISM excitation. In section 6 one- 
phase radiative transfer modeling of the CO line ratios is used to 
systematically extract the average densities, temperatures, and dy- 
namical states of the average ISM environments encountered in our 
sample, and assemble a comprehensive picture of their range. In 
section 7 we present the CO SLED excitation range possible from 
J=l— up to J=7-6, and discuss their power sources. There we also 
discuss the (CO SLED)o(ISM state) degeneracies, and their im- 
pact on the classification of SF modes in galaxies, using the CO 
line ratios of recently discovered near-IR selected gas-rich disks at 
high redshifts as an example. We present our conclusions in Section 
8, where we also list a few important questions regarding the ex- 
pected effects of ISM excitation conditions on molecular gas mass 
estimates in LIRGs and ULIRGs, the subject of our second paper 
(Papadopoulos et al. 2012). Throughout this work we adopt a flat A- 
dominated cosmology with Ho=71 kms -1 Mpc -1 and f2 m =0.27. 



2 THE SAMPLE 

The sample was drawn from two CO J=l-0 surveys of LIRGs by 
Sanders et al. 1991, and Solomon et al. 1997 (themselves drawn 
from the IRAS BGS flux-limited sample with feo M m>5.24 Jy; 
Soifer et al. 1987, 1989; Sanders et al. 2003), so that all galaxies in 
our sample have at least one measurement of this basic H 2 -tracing 
line. We imposed two additional criteria, namely: a) z ^ 0.1 (the 
maximum redshift for which the JCMT B-band receivers can tune 
to CO J=3-2), and b) compact CO-emitting regions (sizes from 
CO interferometric images) so that single telescope pointings or 

The James Clerk Maxwell Telescope is operated by the Joint Astronomy 
Centre (JAC) on behalf of the Science and Technology Facilities Council 
(STFC) of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada 
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Figure 1. The redshift distribution of the combined sample (see section 2). 



small maps can record total line fluxes. If CO images were unavail- 
able then radio continuum (Condon et al. 1990, 1996; Crawford et 
al. 1996), sub-mm dust emission (Lisenfeld, Isaak, & Hills 2000; 
Mortier et al. 2011) or near-IR images (Zenner & Lenzen 1993; 
Murphy et al. 1996; Scoville et al. 2000) helped place an upper limit 
on the size of their CO-bright emission unaffected by extinction. 
This is because sub-mm emission, sensitive to both the warm dust 
associated with molecular gas and the cold dust concomitant with 
the more extended HI (Thomas et al. 2001) sets su b-mm>0co, 
while the radio continuum, tied to the far-IR emission through a 
well-known correlation, traces the star forming regions of LIRGs 
(Condon et al. 1990) where most of the molecular gas lies. 

We have also extensively searched the literature for all local 
LIRGs whose total CO J+l — >J fluxes from J+l=l up to J+l=3 are 
measured (for higher-J lines such measurements are currently non- 
existent in the local Universe). Unfortunately this forced the omis- 
sion of many, often beam-matched, observations (e.g. Devereux et 
al. 1994; Mauersberger et al. 1999; Mao et al. 2011) as in many 
such cases CO line emission extends well beyond the nuclear re- 
gion where such observations were conducted. The combined sam- 
ple of 70 LIRGs (Table 3) is currently the largest for which their 
molecular gas can be studied using CO lines which can probe phys- 
ical conditions from the quiescent to the star forming phase. Its 
redshift distribution is shown in Figure 1, while its IR luminos- 
ity (the SF-powered part) range is L^~(10 10 -2x 10 12 )L Q (com- 
puted over A=(8-1000) (im). The sample is sparce towards low IR 
luminosities since galaxies with Lir~(1-5)x 10 10 L Q are usually 
extended (e.g. Tinney et al. 1990) and thus rarely have total CO 
J+l— >J fluxes available for J+l ^3. We tried to alleviate this bias 
by including the few Iow-Lir galaxies whose total CO J=l-0, 3-2 
line luminosities have been recently measured (Leech et al. 2010). 



3 THE OBSERVATIONS 

We used the A 3 receiver (21 1-276 GHz, DSB operation) on the 15- 
meter JCMT on Mauna Kea in Hawaii (USA), to observe the CO 
J=2-l (230.538 GHz), and 13 CO J=2-l (220.398 GHz) lines, and 
its B3 (315-370GHz) and W/C (430-510GHz) receivers operat- 
ing single sideband (SSB) for CO J=3-2 (345.796 GHz) and J=4-3 



(461.041 GHz) line observations in our sample. The observations 
were conducted during several periods from 1999 up to 2010 (see 
Table 1 for specific periods and typical system temperatures). CO 
J=3-2 observations beyond 2008 utilized the new 16-beam HARP- 
B SSB receiver (325-375 GHz). The decomissioning of the W/C 
JCMT receiver before completion of the survey as well as the sev- 
eral CO J=4-3 lines redshifted into the deep 450 GHz atmospheric 
absorption band meant that such measurements could be conducted 
for only 10 out of the original sample of 30 LIRGs. 

The Digital Autocorrelation Spectrometer (DAS) was used 
in all JCMT observations until 2006, while the new spectrom- 
eter ACSIS was employed aftewards. At 345 GHz we used 
its 920MHz (~800kms _1 ) or 1.8GHz (~1565 kms" 1 ) band- 
width mode, depending on the expected line width and the 
need for maximum sensitivity (i.e. when 920 MHz bandwidth 
was sufficient to cover the line, dual-channel operation was pos- 
sible with B3, and was used for better sensitivity). For the 
high frequency W/C observations the widest 1.8 GHz band- 
width was used throughout whose ~1170kms _1 velocity cover- 
age adequately covered the FWZI of all the CO 4-3 lines ob- 
served. For the CO, 13 CO J=2-l lines both bandwidth modes 
were used, yielding ~( 1200-2345) km s _1 velocity coverage. 
Beam-switching with frequencies of ^ c ho P =(l-2) Hz, at throws 
of 120"-180" (in azimuth) ensured flat baselines. The beam 
sizes were: e HP Bw(230 GHz)=22", 9hpbw(345 GHz)=14" 
and 9hpbw(461 GHz)=ll". We checked and updated the point- 
ing model offsets every hour using continuum and spectral line 
observations of strong s ources, with average residual pointing rms 
scatter a T =^a% + aL<2.5". 

3.1 The CO J=6-5 observations 

The first measurements of the CO J=6-5 line (691.473 GHz) were 
conducted for the luminous ULIRG/QSO Mrk231 and the LIRG 
Arpl93 in our sample using the old JCMT W/D band (620- 
710 GHz) receiver (operating in SSB mode) on 20 of February 
and 22 of April 2005 respectively, under excellent, dry condi- 
tions (T220GHz<0.035). The typical system temperatures were 
T sys ~(3700-5500) K (including atmospheric absorption). The 
DAS spectrometer was used in its widest mode of 1.8 GHz (~ 
780 km s" 1 at 690 GHz), and beam switching at frequencies of 
^cho P =2 Hz with azimuthal throws of 60" resulted in flat baselines. 
The beam size at 691 GHz was 6hpbw=8". Good pointing with 
such narrow beams is crucial and was checked every (45-60) mins 
using differential pointing with the B3 receiver (350 GHz). This al- 
lows access to many more suitable compact sources in the sky than 
direct pointing with the W/D receiver at 690 GHz, and was found 
accurate to within cr r ~2.6" (rms) during that observing period. 

The other CO J=6-5 measurements were conducted during 
2009, with the upgraded W/D receiver equipped with new SIS mix- 
ers (effectively the same type installed at the ALMA telescopes 
in this waveband) which dramatically enhanced its performance. 
The resulting low receiver temperatures (T rx ~550 K) allowed very 
sensitive observations with typical T sys ~(1500-3000) K (includ- 
ing atmospheric absorption) for T220 GHz~0.035-0.06. Dual chan- 
nel operation (after the two polarization channels were aligned 
to within <1") further enhanced the W/D band observing capa- 
bilities at the JCMT (see Table 1). The ACSIS spectrometer at 
is widest mode of 1.8 GHz was used, while in a few cases two 
separate tunings were used to create an effective bandwidth of 
~3.2GHz (~1390kms _1 at 690 GHz) so that it adequately cov- 
ers (U)LIRG CO lines with FWZI~(800-950) km s" 1 . Rapid beam 
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Table 1. JCMT: observing periods, CO lines, typical system temperatures 



Year 


periods 


spectral lines 


T S ys (K)" 


1QQQ 

1777 


06/07 90/07 


rn 3 9 A 3 


4S0 700 (TK X\ 1 1 ^0 9^80 fW/Ct 


2001 


10/1? 98/1? 

1 u/ 1 z. — Z.O/ 1 z. 


CO 3—2 


^80-SSO fR ^ 


2002 


94/09 9^/09 


CO 3—2 


6^0-6Q0 (Ti 1} 




1 7/04—18/04 


CO 3—2 


760 1 100-1700 fR ^ 




17/06 30/flfi 


CO 9 13 2 

Z. — 1, J — _ 


SSOfA^ 4^0-6^0 Q00-1900fR^ 




90/1 1 9V1 1 


CO 4 3 


9600-4800 6000-Q^OO (W/C\ b 


2003 


06/1 1 


CO 3—2 


4000-4400 (TK T\ 


2004 


20/01 


CO 3—2 


^0-690 (A X\ 




02/04— 29/0*5 


12,13 pfl 91 pQ o_ 7 

V V/ Z, 1 . V V / . ' Z. 


300-450 I A 31 470-2400 CR 31 6 


t( 


1^/07 9S/08 


12,13pn 9 1 rn q 7 
V V 7 z. — 1, V V / J — z. 


^90-S10 CA 800-1000 fR ^"1 




98/OQ 10/1 1 


I2,i3pn ? 1 m 3 ? 

V . V 7 Z. — 1, V V / J — Z. 


^00-600 f A ^0 fR ^ 

JUU-UUU \r\ J J, JjU- 1 JJU \D D ) 


(t 


18/11 


CO -4-3 


94S0-9800 fW/O 


2005 


20/02 


CO 6-5 


4500-5200 rw/D^ 




17/04-9^/04 


ro ^ 9 ^ s 

V. V 7 J — Z., *t J, U — J 


790-1 1 00 fR ^"1 1 600 (WIC\ ^700-SSOO ^W/D"! 

/ Z.U 1 1 L/VJ \L> -J )i 1 UVJVJ ^ VV / »— 7' ■-* * W -/JW V VV / ^ 




22/08 


13 CO 2-1 


41 5-470 (A 3"! 


« 


10/10-28/10 


12 - 13 CO 2-1 


230-420 (A 3) 




15/12-31/12 


12 - 13 CO 2-1, CO 3-2 


330-490 (A 3), 750-980 (B 3) 


2006 


15/12-18/12 


13 CO 2-1 


250-280 (A 3/ACSIS) 


2007 


16/12-21/12 


13 CO 2-1 


230-280 (A 3/ACSIS) 




22/02-24/02 


13 CO 2-1 


310-320 (A 3/ACSIS) 


2008 


09/05 


CO 3-2 


1630 (HARP-B/ACSIS) 


2009 


06/01-07/01 


CO 6-5 


2600-3200, 7500 (W/D, ACSIS) 6 




22/01-25/01 


CO 6-5 


1700-2600, 9000 (W/D, ACSIS) 6 




01/02 


CO 6-5 


2200-3500 (W/D, ACSIS) 




27/01 


CO 6-5 


2000-2900 (W/D, ACSIS) 




02/03 


CO 6-5 


1300-1500 (W/D, ACSIS) 




13/03-15/03 


CO 6-5 


1900-3100 (W/D, ACSIS) 


2010 


11/09 


CO 6-5 


1400-1800 (W/D, ACSIS) 



a T sys values include atmospheric absorption. High values of T sys (B 3)>600 K and 
T sys (A3)>350K) are measured in a few cases close to tuning range limits (K 5 ] 5y <320GHz for 
B 3, Kjky ^215 GHz for A 3), or when \ u s ^ y — 325 GHz|<5 GHz, i.e. close to a strong atmospheric 
absorption feature at 325 GHz. 

b The high T sys values found only for VIIZw 31, a circumpolar source at the JCMT latitude, 
observed at elevation of ~ 30° 

Table 2. Point source /T (Jy/K) conversion factors and HPBW values used for the observations and literature data 



Telescope" 110-115GHZ 6 210-230GHz 315-345GHZ 6 430-461 GHz 6 620-710GHz 6 



IRAM30-m (S„/T^) 


6.3 C (22") 


7.9, 8.7(11") 








JCMT(S„/Tp 




25-28(22") 


28-38(14") 


50-74(11") 


49-62(8") 


NRAO 12-m(S^/T* ) 


35 d (55") 


55(32") 








FCRAO 14-m (S„/T^) 


42(45") e 










Onsala20-m (S„/T^) 


31(33") 










SEST (S„/T mb ) 


27(45") / 


41(22"/ 








HHSMT (Si//T m t,) 






50 9 (23") 






NR0 45-m (S„/T mb ) 


2.4 h (14.5") 










CSO(S,/T mb ) 






43 l (21") 







"The telescope and temperature scale type (see Kutner & Ulich 1981 for definitions). 
b The frequency range per receiver. 

c For T^=(B eff /F off ) X T mb =0.789 X T mb : S„/T mb =4.95 Jy/K, used to obtain the CO 1-0 fluxes from all the 
IRAM 30-m spectra in the literature that are reported in the T mb scale (unless a different S„/T mb is mentioned). 
d Measured at 1 10 GHz (section 2.3.2), also in the NRAO 12-m User's Manual 1990 Edition, Figure 14. 
e The Tj^ scale is sometimes used to report data from the FCRAO 14-m (e.g. Sanders et al. 1986), for which 
S„/T^ = rite (S„/TX) = 31.5(Jy/K) (for r, fss = 0.75) is adopted. 

f http://www.ls.eso.org/lasilla/Telescopes/SEST/html/telescopeinstruments/telescope/index.html 
9 The Heinrich Hertz Submillimeter Telescope (Arizona, USA) (from Narayanan et al. 2005). 

h For the NRO 45-m telescope in Nobeyama (Japan) at 115 GHz: V = (8k B /7rD 2 )(r? mb /r; a ) and adopting 
rj a = 0.32, r] mh = 0.44 (from the NRO website). 

^altech Submillimeter Observatory (CSO): V = (8k B /7rD 2 )(»7 mb /»y a ), where D = 10.4m, r\ mb = 0.746 and 
01/2 = 1.22(A/D) yielding Tja/Vmb = 0.76 (using Q A A e = A 2 and fj a = A c /A g , r) mb = Q mb /n A - 
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switching at i^ c ho P =4Hz (continuum mode) and azimuthal throw 
of 30" yielded very flat baselines under most circumstances. The 
pointing model was updated every 45-60 mins using observations 
of compact sources with the W/D receiver, as well as differential 
pointing with the A 3 receiver, yielding rms residual error radius of 
cr r ~2.2". The final CO J=6-5 observations were conducted in 201 1 
during which I Zw 1 was observed, with only one W(D) receiver 
channel functioning, under dry conditions (Y220 GHz<0.05) that 
yielded T sys ~( 1400-1 800) K. The same beam-switching scheme 
was used, while two separate tunings yielded an effective band- 
width of ~3.2 GHz covering the wide CO line of this ULIRG/QSO 
(e.g. Barvainis et al. 1989). The pointing uncertainty remained 
within the range of previous CO J=6-5 observations. Nevertheless 
we wish to note that isolated cases of large pointing offsets reduc- 
ing the observed CO J=6-5 line fluxes have been found (e.g. for 
Arp 220, see PlOa and PlOb), and may have affected a few of these 
highly demanding CO line measurements. 



adopted the JCMT measurement (as its wider beam is less prone to 
flux loss due to pointing offsets and/or beam-throw/flux-loss uncer- 
tainties) if a discrepancy larger than the aforementioned was found. 



4 DATA REDUCTION, LINE INTENSITY ESTIMATES, 
LITERATURE DATA 

In both telescopes the output spectra are in the scale (see Kut- 
ner, & Ulich 1981). We inspected all individual 10min/(4-6)min 
JCMT/IRAM spectra for baseline ripples and to clip any intensity 
"spikes" in individual channels. The edited spectra were then co- 
added using a 1/ct 2 -weighting scheme and linear baselines were 
subtracted from each final co-added spectrum. These spectra are 
shown in Figures 2 and 3 and were used to derive the velocity- 
integrated molecular line flux densities from 



3.2 The IRAM 30-m observations 

Observations of CO, 13 CO J=l-0, 2-1 with the IRAM 30-m 
telescope were conducted during two sessions in 2006 namely, 
June from 20 to 25, and November from 26 to 28. In both pe- 
riods the A100/B100 (3 mm) and A230/B230 (1mm) receivers 
were used, connected to the 1 MHz (A100/B100, 512 MHz) and 
4 MHz (A230/B230, 1 GHz) filterbanks. During the first period the 
A230/B230 receivers were used to observe the 12 CO J=2-l line. 
If the latter was strong and detected in about a hour or less, the 
1mm receivers were then re-tuned to 13 CO J=2-l. For sources 
with very weak 12 CO lines (e.g. 08030+5243, 08572+3915), no 
attempt of observing 13 CO J=2— 1 was made. Data were acquired 
under New Control System (NCS) in series of four-minute scans, 
each comprised of eight 30-sec subscans. The typical system tem- 
peratures (including atmospheric absorption) for the CO 2-1 ob- 
servations were T sys (210-230 GHz)~(220-500) K, with the lowest 
mostly during the 13 CO observations (though for occasional tun- 
ings towards the edge of the band and/or bad weather conditions 
T sys ~(700-900) K). For most sources data were acquired in two or 
more different days to ensure a line detection, and as a consistency 
check. Pointing and focus were checked frequently during the ob- 
servations with residual pointing errors <j r ~3" (rms). 

During the November period receivers A100/B100 were used 
to observe 13 CO J=l-0 line simultaneously to the 12 CO J=2-l 
line observed with A230/B230 (tuned to the same line each time). 
The pointing error stayed <3" (rms), except during November 26 
when it went up to ~6" (corresponding data were omitted). The 
typical system temperatures were T sys (110 GHz)~l 10-160 K, 
T sys (115GHz)~(200-380)K, and T sys (210-230GHz)~(330- 
425) K. Finally, in order to maintain very flat baselines, the 
wobbler switching (nutating subreflector) observing mode with 
a frequency of 0.5 Hz and beam throws of 180"-240 was 
employed during both observing sessions. The beam sizes 
were: H pbw(11O GHz)=22" and e H PBw(210-230GHz)=ll", 
with corresponding beam efficiencies^ of: B e ff(110 GHz)=0.75, 
B eff (230 GHz)=0.52 and B cft (210 GHz)=0.57, and forward beam 
efficiencies of F e ff(3 mm)=0.95 and F c fr(l mm)=0.91. We also 
note that in most cases we had redudant CO J=2— 1 measure- 
ments with the JCMT, and then: a) adopted the average when 
JCMT/IRAM values agreed to within 20% (most cases), or b) 

2 http://www.iram.fr/IRAMES/index.htm 



Siinc = / S v dV = - s JL-K B {x) / TXdV 

J AV Va^L) J Av 

= V -^K c(x) { T XdV , (1) 

r > a J AV 

where rjcMT=15.62, riRAM=3.905 and rf a is the aperture effi- 
ciency defined against the scale (ril=i]a/Vrss, where r\ a is the 
aperture efficiency measured against the T^ scale, as is more typ- 
ical, and r/ rss is the rearward spillover and scattering efficiency, 
Rohlfs & Wilson (1996), Eqs 8.16, 8.17). The factor K c (x) = 
x 2 /(l — e _x ), with x=8 s /(1.28hpbw) and # s =source diame- 
ter, accounts for the geometric coupling of the beam (its gaussian 
part) to a disk-like source, when a CO emission size was avail- 
able, and 2a r <8 s <0uPBW where oy is the pointing error radius 
(see 4.1). The total point-source conversion factors S„/T adopted 
for the JCMT, the IRAM 30-m telescope, and all the data gleaned 
from the literature (for the corresponding output antenna tempera- 
ture scales) are comprehensively tabulated in Table 2. 



4.1 Aperture efficiencies, line intensity uncertainties, biases 

Aperture efficiencies of large high-frequency sub-mm telescopes 
such as the JCMT can change significantly (especially for 
i^>460 GHz) depending on a variety of factors (e.g., elevation, ther- 
mal relaxation of the dish or its re-shaping after an holography ses- 
sion). In order to track them over an decade of observations (during 
which the JCMT dish has been re-adjusted quite a few times) we 
conducted frequent aperture efficiency measurements using planets 
and adopted the average 77* obtained per observing period for de- 
riving the line fluxes of all the sources observed during it. In many 
cases, as a cross-check, we distributed the measurements of very 
CO-luminous LIRGs over several widely separated periods, during 
which very different aperture efficiencies (sometimes up to a factor 
of ~2) were often measured. In all such cases Equation 1, with the 
appropriate 77* values, yielded velocity-integrated line fluxes in ex- 
cellent agreement. Indicatively most aperture efficiencies measured 
for the JCMT lay within ?7*~0.41-0.56 (B-band, 315-350 GHz), 
and Va~0.21-0.31 (C-band, 430-461 GHz). For the three periods 
of the more demanding W/D band observations we derived 7^=0.25 
(2005), 0.32 (2009), 0.27 (2010) from planetary measurements. 
The uncertainties for the reported velocity-integrated line flux den- 
sities have been computed from 
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Figure 2. The high-J CO J+l— »J, J+1^3 spectra. In the few cases where all three CO J=3-2, 4—3 and 6-5 lines are available we omit the overlay of CO 
J=2-l in order to reduce confusion (the J=2— 1 lines are all shown in Figure 3). The velocities are with respect to V pt=cz C o(LSR) (Table 3), and with typical 
resolutions AV c v 1 ~(10-50)kms _1 . A common color designated per transition is used in all frames. 



Slinc 



V -Lav /th V Iav /cai \ V J 



-, 1/2 



(2) 



where <5Tav is the stochastic error of the average line intensity 
Tav (averaged over the line FWZI AV), r\ is the telescope effi- 
ciency factor (used to derive the integrated line flux from the tem- 
perature scale of the output spectrum, e.g., rf a for T^IE and &r\ its 



uncertainty. The first term is estimated from the spectra shown in 
Figures 2, 3 using 



V Tav / th Tav V NAvNbas / 



(3) 



where 5T c h an is the stochastic intensity dispersion, estimated 
from the line-free part of the spectrum (for a given ve- 
locity channel width AV c han), while NAv=AV/AV c han and 



3 More precisely the 3nd term in Equation 2 embodies the uncertainty 
of total efficiency factor TjiTjf sa needed to convert the scale (cor- 



rected only for atmospheric extinction) to a final fully corrected (apart from 
source-beam coupling) (or T m j,) scale. 
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Figure 2 - continued The high-J CO J+l— >J, J+1^3 spectra. In the few cases where all three CO J=3-2, 4-3 and 6-5 lines are available we omit the overlay 
of CO J=2-l in order to reduce confusion (the J=2-l lines are all shown in Figure 3). The velocities are with respect to V pt=cz c0 (LSR) (Table 3), and with 
typical resolutions AV c h~(10-50)kms" 1 . A common color designated per transition is used in all frames. 



Nbaa=2AVbaa/AV c han are the number of channels within the line 
FWZI and the line-free baseline (with AVbas/AV c han channels 
symmetrically around the line) respectively. The second term in 
Equation 2 accounts for line calibration errors (due to a host of fac- 
tors such as imprecise knowledge of the calibration loads, uncer- 
tainties in the atmospheric model and the derived extinction etc). 
Observations of numerous strong spectral line standards and plan- 
ets during each observing period yielded intensity dispersions of 
~15% (230 GHz and 345 GHz), ~20% (460 GHz), and ~25% 
(690 GHz), which we adopt as the combined calibration (cal) and 



Sri/rj uncertainties per observing band at the JCMT. For the 30-m 
we consider these to be ~15% for both 3 mm and 1 mm bands. 

Finally, even with the accurate tracking and pointing achiev- 
able by enclosed telescopes such as the JCMT, the residual rms 
pointing errors and the narrow beams of large mm/sub-mm tele- 
scopes at high frequencies can lead to a substantial and systematic 
reduction of measured fluxes of compact sources. We try to account 
for this as described in PlOa, by applying a (G) scaling factor to 
the line fluxes of all point-like sources with CO emission region di- 
ameters of 6? co <2a r (where cr r (JCMT)=2.5" and cr r (IRAM)=3" 
are the pointing error radii). This factor is (see PlOa): 
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Figure 2 - continued The high-J CO J+l— >J, J+1^3 spectra. In the few cases where all three CO J=3-2, 4-3 and 6-5 lines are available we omit the overlay 
of CO J=2-l in order to reduce confusion (the J=2-l lines are all shown in Figure 3). In the case of NGC 5256 only the CO J=l-0, 3-2 lines are available. The 
velocities are with respect to V op t=cz co (LSR) (Table 3), and with typical resolutions AV c h~(10-50) km s — 1 . A common color designated per transition is 
used in all frames. 



G(<7r) = 1 



!ln2 



/2 



(4) 



where 9y 2 is the beam HPBW, ay / y/2 is the rms pointing error per 
pointing coordinate. For 9 co ^2a T , K c (x) is replaced in Equation 
1 by G(cr r ) as the beam-source coupling correction is overtaken 
by the pointing error correction. At 345 GHz and 460 GHz for the 
JCMT we obtain (G 345 ) = 1.087 and (G 46 o)=l-15, (with neglible 
correction at 230GHz), while for the 30-m (G23o) = l-20 (and neg- 
ligible correction at 115 GHz). The CO J=6-5 observations with 
HPBW~8" are those most succeptible to this bias and (Gs9o) has 



been estimated from the pointing rms per observing session, yield- 
ing a range (G69o)=l- 17-1.37. For sources with CO (or sub-mm 
dust emission) sizes of 2a r <9 s <8 1 / 2 the K c factor is used in Equa- 
tion 1 . Finally in the cases where large offsets were found between 
the presumed CO source center and the observed positions in the 
literature (see discussion in 4.2) we applied a beam-shift correction 
factor of K s h=exp [4 ln2 (A#/0hpbw) 2 ] where AO is the (beam 
center)-source offset (see Table 4). 
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Figure 2 - continued The high-J CO J+l— >J, J+1^3 spectra. In the few cases where all three CO J=3-2, 4-3 and 6-5 lines are available we omit the overlay 
of CO J=2-l in order to reduce confusion (the J=2-l lines are all shown in Figure 3). The velocities are with respect to V pt=cz c0 (LSR) (Table 3), and with 
typical resolutions AV c h~(10-50)kms^ 1 . A common color designated per transition is used in all frames. 



4.2 Incorporating data from the literature, the final dataset 

A detailed literature search for all total CO and 13 CO line fluxes 
available for LIRGs enlarged our sample (see Table 3), while allow- 
ing also a consistency check using the duplicate measurements per 
object (especially for the 1=1-0 transition). In most cases we find 
good agreement among the various CO J=l-0 fluxes reported in the 
literature, and between the (much fewer) reported J=2-l, 3-2 line 
fluxes and our measurements, within the expected uncertaintiefl 



Most cases of serious discrepancies among CO line fluxes were 
rectified after accounting for the different source positions often 
used. This typically occured because of the positional uncertain- 
ties inherent in the optical identification of heavily dust-obscured 
sources within the large IRAS position error ellipse (e.g. Solomon 
et al. 1997) and/or the better CO positions obtained for some LIRGs 
(e.g. via interferometry) and used for subsequent CO observations 
after the original J=l-0 detections with usually large (~45"-50") 



4 In the literature the uncertainties of mm/sub-mm line measurements are 
often underestimated, with only the thermal rms error reported (the first 



term in the expression in Equation 2). In all such cases we assumed a 15% 
of calibration and r\ a uncertainty in addition to the one reported. 
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Figure 3. The CO, 13 CO line data. The velocities are with respect to Vopt=cz co (LSR) (Table 3), and with typical resolutions AV c h~(35-90) km s 1 . A 
common color designated per transition is used in all frames. 



beams were made (e.g. Sanders et al. 1991; Young et al. 1995). 
Moreover, even when a particular LIRG is optically identified, the 
multi-component/interacting nature of many such systems can con- 
found the choice of a pointing center for CO single dish observa- 
tions (e.g. Leech et al. 2010). In the cases where we found different 
positions used for CO observations reported in the literature, we 
"shift/scale" the corresponding line fluxes to a common position, 
the interferometrically-derived CO emission peak of the LIRG (or 
its near-IR and/or radio continuum peak if the former was unavail- 
able) assuming a source much smaller than the beams used (the 
case for most LIRGs). This rectified many of the discrepant CO 
line fluxes, with only a few inconsistent ones remaining, mostly be- 



tween new data and much older CO observations (mainly those re- 
ported in the Sanders et al. 1991 CO J=l— survey). These are likely 
the result of improved calibration techniques (e.g. good sideband 
rejection), pointing accuracy, and faster beam-switching available 
for mm/sub-mm telescopes now than in the past. 



4.2.7 Resolving a large discrepancy: new CO J =1-0 
observations of IRAS 05189-2524 

The largest remaining CO line flux discrepancy is for the LIRG 
IRAS 05189-2524 with Si =(91±18) Jykms- 1 obtained for CO 
J=l with the NRAO 12-m telescope (Sanders et al. 1991) and 
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Figure 3 - continued The CO, 13 CO line data. The velocities are with respect to Vopt=cz C o(LSR) (Table 3), and with typical resolutions AV c h~(35- 
90) km s - 1 . A common color designated per transition is used in all frames. 



Sio=(52 ± 10) Jykms -1 with SEST (Strong et al. 2004). The 
larger value would place this ULIRG to the lowest end of gas exci- 
tation of the entire sample, with r32 =0.24 and r2i=0.32, typical for 
the coldest most quiescent GMCs found in the Galaxy and M31 
(Loinard et al. 1995; Allen et al. 1995; Fixsen et al. 1999). We 
used the 12-meter telescop^fl to re-observe the CO J=l— line in 
this ULIRG on the nights of 12, 13 and 24 March, 2008. The fi- 
nal co-added CO J=l-0 spectrum (Figure 4) was used to estimate 
a velocity-integrated line flux of Sio=(43 ± 8) Jy km s _1 , in good 



agreement with the value reported by Strong et al. 2004 but less 
than half that reported by Sanders et al. 1991. We report the aver- 
age of ours and the Strong et al. 2004 value (Table 4). 



5 The Kitt Peak 12 m telescope is operated the Arizona Radio Observatory 
(ARO), Steward Observatory, University of Arizona 
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Figure 3 - continued The CO, 13 CO line data. The velocities are with respect to V pt=cz c0 (LSR) (Table 3), and with typical resolutions AV c h~(35- 
90) km s~ 1 . A common color designated per transition is used in all frames. 



5 THE MOLECULAR GAS IN LIRGS 



The ratios of velocity/area-averaged brightness temperatures of 
CO lines provide the excitation indicators of the average state of 
the molecular gas, allowing comparisons among LIRGs and well- 
studied local ISM environments in the Galaxy where individual ex- 
citation mechanisms (e.g. SNR shocks, far-UV/IR photons, cosmic 
rays) can be more easily identified. The luminosity of a line (x) 
used to compute such brightness temperature line ratios is 



T b - da dV 



AV J A s 

2 



1 + Z 



S„dV, 



(5) 



where Tb, x is the rest-frame brightness temperature of the line, 
AV, A a are the line FWZI and source area, Dl the luminosity dis- 
tance, and S„ the line flux density. After substituting astrophysical 
units 



3.25XHT 



D£(Mpc) 



1 + z 
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100 GHz J 



Jy kms -1 



,(6) 



where L x is in Kkm s pc and v Xirea t ls the rest frame frequency 
of line (x). This expression is used to compute the CO line ratios 
with r.j+i j=L J+1 j/L 10 . The conversion to ordinary luminosity 
units (L©), used to express the total line luminosities (L x = J L„ dv) 
in CO SLEDs is 



ij X — T, 



3.18 x 10 



4 ( 



,rest 

100 GHz, 



r 






10 9 L, 



(7) 



Kkms : pc 2 . 



where Li 

In Tables 7 and 8 we give the CO line ratios for our sam- 
ple, and in Figures 5, 6, 7 their frequency distributions. These re- 
veal well-excited lines ((r2i)=0.91, (r32)=0.67) but also a signifi- 
cant excitation range, from low (r2i~0.6-0.7, r32~0.3) to a high- 
excitation phase (i~2i, 32-^1)- Luminous CO J=4—3 lines clearly 
mark the emergence of a second highly excitated phase with 
r 43<;r32 and r32<r2i<l in several LIRGs (Table 7). This could not 
happen if the average state was dominated by one phase with opti- 
cally thick CO lines since then 1 > r/21 > r32 > r43 > ... (higher- J 
lines become fainter as subthermal excitation progressively sets in 
and/or as Ej+i, j/QtsTk) becomes >1). A high-excitation phase 
overtaking the lower-excitation one at the transition where the line 
ratios of the latter fall below unity can indeed reverse this inequal- 
ity series, provided that this new phase has well-excited and opti- 
cally thin CO SLEDs. In some cases this occurs already from J=3-2 
withr 3 2>r 2 i andr 2 i<0.6-l (e.g. IRAS 10190+1322, Arp 220). We 
note that while the emergence of a second highly-excitated phase 
in global CO and 13 CO line emission has been known for some 
time (e.g. Aalto et al. 1995 and references therein), indications for 
well-excited and partially optically thin global CO SLEDs are rare. 
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Table 3. Observational parameters of the combined LIRG sample 
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0.0445 


(194.5) 




in" 


o"1 
i u ) 


R1 " (m\ 


jp,m 


00177 0840* rMnr 1 S7i 


00 14 46 48 


08 21 47 8 


0.0055 


(23.3) 




("0" 

( u 


0"1 

u ) 


OvJ ^Lvj,SII1,a J 


2P,m 




00 51 14 92 

\J\) JT. 7i- 


-1-17 41 15 5 

1^ 1 Z. tl JJ.J 


0.0611 


(270.3) 




("0" 


0"1 
u ) 


< i2 // (co) 


3 m 4 p 


01 0^1 174^* ("Arn916^ 
i./ 1 ujj — i / t^u v*" M z. jvj 7 


01 07 47 00 


— 1 7 10 74 

— i / jvj z.^t.vy 


0.0200 


(85.8) 




(n" 
v u 


0"1 
u ) 




5 P 6 m 


01077 1707* 

VJ1 VJ / / — 1 1 \J I 


01 10 08 70 


— 16 S1 11 

1UJ1 11 .V/ 


0.0351 


(152.3) 




(0" 


0"1 


< 1 'V'f'srrO 

-LiJ ^Slliy 


5 p 6 m 


014184-16S1* fTTT7wl*i"l 


01 44 ^0 SO 


-1- 1 7 06 OS 
i / vjvj vjo.vj 


0.0274 


(118.2) 




(0" 


ft"\ 


< 1 V'Csm'l 

-LiJ ^Slliy 


5 P 6 m 


0207 1+3857 fNGC 828 VI 7w 1 771 


07 1 09 4^ 


-1-19 1 1 26 3 

| J 1/ 1 1 Z. V7 . J 


0.0178 


(76.2) 




(0" 


ft"\ 


14"(co) 


7P,m gm 


07080+3775* fNGr8341 

UZLKM/+ J / Z J ^1>IVJ\^. O J^ J 


02 1 1 01 55 


-1-17 40 01 1 


0.0154 


(65.8) 


( 


-2.5" 


t —2.0") 


lO ^uiiiy 


9 P 10™ 


0211 4+04 56* fMrk 10271 

VJZ. I 1 ^ l V J ^ J V7 I IV 11 IV 1 VJZ. / J 


02 1 4 05 60 


-Lf) ] 1 77 7 
n^vj j i vj £.1,1 


0.0297 


(128.3) 




(0" 


0"1 


1 7"(sm xl 

_L i v^aiii, a ) 


5 p 6 m 


02321 OQOfl fNfiT 985 Mrlc 10481 


02 14 17 74 


—08 47 14.7 


0.0430 


(187.7) 




(0" 


0"1 






02401 0013* ("NGC 10681 

UZ+Ol — U\J 1 J ^INUV iouo_j 


07 47 40 74 

V/Z. T'Z. T^l.J. / T 


—00 00 47 6 

UV' V'U t / .u 


0.0037 


(13.3) 




(0" 


0"1 
u ) 


/dO" i nc\ cm y \ 
L ±v) I v Uvj,M11,A ) 




02483+4302 

VJZ.T-0 JTH J VJZ. 


07 51 16 01 

V/Z. J 1 JU.U 1 


-1-41 1 S 10 8 
n^^+j ij iv^.o 


0.0514 


(225.8) 




(0" 


0"1 
u ) 


y . I o ) 


jp,m 


07S174-1446* runr 1 7i6Q't 

UZ. J 1 Z.T 1 t-t-VJ ^ UUv, iJU7j 


07 S4 01 80 


J-14 S 8 14 

It JO li.U 


0.0312 


(135.0) 




(0" 


0"1 
u ) 


i -^0"/ cm r*m v ^ 
OU ^6111,^/111, A ) 


5P (^m ^ m 


011SQ4-1 S71 


0^ ^4n Q0 


-4-1 S 19 SS 
i J J z. j j . vj 


0.0353 


(153.2) 




(0" 


0"1 


A WrrrO 
'-t.O V *-lll / 


5 P 13 m 


047174-1416 


04 76 04 Q4 
yj^t z.u ui.77 


+ 14 43 37 9 


0.0796 


(356.4) 


(- 


-1.13' 


+0 4"1 


ft " ( r*m c m v ^ 


14 p 6 m 13 m 


05083+7936 fVll 7w 031 1 


05 1 6 46 5 1 


_I_7Q 40 1 7 S 

/ y taj i^.j 


0.0543 


(239.0) 




(0" 


0") 


1 3"(co1 

_ . ' ) V UVJ J 


jp,m 


05189-2524 

UJ 107 Z — IZ.'-T 


05 2101 11 

v/ j — y vi.ii 


—25 21 45.9 


0.0427 


(186.4) 


(+4.0" 


4-1 0"1 


< "\"(cm irl 


9 P 15 m 13 m 


08030+5243 

VJOVJ J 


08 06 50 10 


-1-52 35 05 4 

T-'Z JJ VJJ.T^ 


0.0835 


(375.4) 




(0" 


0") 


< 8"fir1 


14 p ,16 m 


08354+2555* (NGC2623 Arn2431 


08 18 24 10 

V/O JO Z/-r.lW 


4-25 45 16 5 

^^Z — ' i J 1V7.J 


0.0185 


(79.3) 




(0" 


0"1 


i.vjo ^v-uy 




08572+391 59 

VJO J 1 iTJ7 U 


09 00 25 41 

\jy vt\j z. j . 1 1 


4-39 03 54 1 

1 J T 7 V/J J7.1 


0.0582 


(256.9) 




(0" 


0") 


< 2.1"(co) 


18 p ' m 


09126+4432 (Am 551 

vjy l i-u t tt^j z. \ yj j j j 


09 1 5 54 90 


+44 19 54.4 


0.0399 


(173.3) 




(0" 


ft"\ 
u > 


1 Q" 1 2"fco cm rihll 

y .zf - iz v^vj,L.iii,viuiy 




0Q1704-6114 nine osioi"i 


09 15 51 51 


-U61 71116 

TUl Z.1 1 1 .VJ 


0.0393 


(171.1) 




(0" 


0"1 


i z< \yy,\s\j} 


14 p 15 m 20 m 


09586+1600* fNGC 30941 

Vjy^JvyVf^yVJVJVJ 11 *i V J V Ju7^ J 


1 01 26 00 

1 VJ VJ 1 £j\J ,VJ\3 


4- 1 5 46 14 


0.0080 


(34.0) 


( 


-2.2" 


— 1 0"1 

±.U ) 


1 1 "fir rml 

i y. \ ii,uiii / 


2l p 22 m 13 m 


10039-3338* TIC 25451 

JVJVJ^J y J JJO V^-^ «J } 


1 06 04 50 

yvi v}\j ut. jv/ 


—33 53 03 

JJ JJ VJJ.V/ 


0.0341 


(147.9) 




(0" 


0") 


< 1 ^"fsml 


5 P 6 m 


10035+4852 

1 VJVJ^J^J 1 1 O J Z. 


1 06 45 81 

1 VJ VJVJ T^J.CjJ 


4-48 37 46 1 

| lu J / T^VJ. 1 


0.0648 


(287.5) 


(+5.2" 


4-2 2"1 


< 1 Vfsml 

-^vj \jyyyj 


14 p 6 m 


1(1173+0878 


1 20 00 1 9 


-4-08 1 3 34 5 

T^O U J+ — ' 


0.0489 


(214.4) 




(0" 


0"1 


< i \" (rr, mil 

*-* V_VvVJ,L/lll^ 


21P,rn 


10190+1322 


10 21 42 60 

± V7 Z. ± T^Z. . VJV7 


j_ 1 1 06 54 4 

i J vjvj jt.t 


0.0765 


(342.2) 




(0" 


0"1 


1 2" 4"frn rihll 

J..ZJ . '-t ( V-.U,L1U1 J 


24 p ' m 


10356+5345* fNGC33101 


1 18 45 90 

1 V7 JO *tJ.7u 


4-51 10 1 1 7 

T^J JV-* 11./ 


0.0033 


(14.0) 




(0" 


0") 
u ) 


^r\"(cn sm xl 

• J VJ V_*-VJ, Sill, Ay 


2P,m 


10565+7448 


1 59 1 8 1 5 


-1-74 17 14 4 

T^Z.^t JZ. J7.7 


0.0428 


(188.2) 




(0" 


0"1 


1.5"(co) 


jp,m 


1 1 191+1900 (PG 1 1 19+1901 

11 1 y IT 1Z.VJVJ yJT v J 111 7T 1 Z.VJ ^ 


11 71 47 1 7 

11 Z.1 t/.iz 


4-1 1 44 1 8 1 

1 1 ii i o . j 


0.0500 


(219.4) 




(0" 


0"1 


< 5"(co) 


25 p ' m 


11931+1456* (7r2810 TTGr 64361 


11 75 45 00 


-I- 1 4 40 16 


0.0341 


(147.9) 




(+1.5 


II n ll\ 
) u ) 


' ( cyjn cmA 
U ^L-lll,slliy 


21 p 13 m 26 m 


1 1 957+5850 ( Am 2991 


1 1 28 12 45 

1 1 Z. O J^.T^J 


4-58 11 45 8 

| JO JJ T^J.O 


0.0103 


(43.8) 




(0" 


0") 


l^'Veo cm xl 

ijtj \k-VJ,LMll,A / 


27 p ' m ,6 m 


1 7 00 1 +02 1 5 * TNfir 404 5 1 

1 i-\j\J 1TIJZ. 1 J ^l"IVJv- J } 


1 7 07 47 ^0 

1 Z. V7Z. ^tZ 


-1-0 1 S8 18 

\^\) L JO JO.VJ 


0.0066 


(28.0) 


( 


-0.9" 


— 1 7"1 


7" (cm) 


21 p 13 m 


121 12+0305 

1^~VJJVJJ 


1 2 1 3 45 77 

X Z. 1 J tj . / / 


4-02 48 39 3 


0.0727 


(324.4) 


( 


4-3.9" 


4-2 1"1 


< 2" (2 9" Vm rihll 

_j Zi IZ*.^J j \ V_. U , LI L/ 1 } 


9 p ,18 m 


12224 0624* 

1 Z.Z.Z.T 1 — LJVJZ.'T 


1 7 25 03 90 


—06 40 53 


0.0263 


(113.4) 




(0" 


0"1 


< 2" (cm) 


21 p 13 m 


12243-0036* (NGC44I81 

1Z,Z,TJ VAIJU ^11 VJV, ii 1 (J / 


1 2 26 54 70 

1 Z. Z.V7 Jt, /V/ 


—00 52 39 

vjvj jz. ^)y.vj 


0.0073 


(31.0) 




(0" 


0") 


< (cm) 


21 p 13 m 


12540+5708 fMrlc2311 

1 Jt^UTJ / VJO ^IVllJS. Z. J 1 ) 


17 56 14 71 

1 Z. JVJ It.Z 1 


4-56 52 25 1 

TJu j^ ^J.l 


0.0422 


(184.1) 




(0" 


0"1 

u ) 


U.OO \\svJ ) 


jp,m 


13001-2339* 


13 02 52 10 

1 J V7Z. J Z. . IV/ 


—23 55 19 

Z.J JJ 17. VJ 


0.0215 


(92.3) 




(0" 


0") 




5 p ,22 m 


11102+1251* (IMGr 50201 

1 J 1 VJZ.T IZiJ 1 ^1 "1 vJv, JVJZ.VJ J 


11 1 7 Q0 

lj lz jy.yv) 


-1-12 35 59 
^i^ jj jy.\) 


0.0112 


(47.7) 


( 


-1.0" 


5 —0.9") 


1 2" (rrr\\ 
i Zi v^uiii j 


21 p 13 m 


Arn718* iTTnr0811^ 

rtljJ Z.JO ^UVJV^. UOJJJ ) 


n 1 s io 70 

lj 1J JI/.Z.V' 


4-62 07 45 


0.0315 


(136.3) 




(0" 


0"1 
u ) 


ell (r>c:ll\( rrn c \ui\ 
o iou n L,iii,u.Ljiy 


5 P 10™ 


1 11 814-1471 TArn 1 Qlt 


11 70 IS 17 

lj Z.VJ JJ.JZ 


4-34 08 77 7 

yj7 VJO Z.Z..Z. 


0.0233 


(100.2) 




(0" 


0"1 


l.S'^co) 


jp,m 


1 i 1 88+0016* rivrir' s 1 04~t 

1 J 1 OOTUV/Ju ^ 1 "1 VJV_ J 1 VJt^ ^ 


1171 71 10 

1 J Z. 1 Z.J . IV 


-1-00 70 17 
n^/vj z.vj jz..vj 


0.0186 


(79.7) 




(0" 


0"1 


9 ^"(rm} 

Z . O V Ulll / 


21 p 13 m 


1 3229-2934 fNGC 51 351 

l j z. z. y £*y *j"~x 11 *< v i v j i jj i 


1 3 25 43 97 

1 J Z< J J . .7 / 


—29 50 01 3 

Z. .1/ J VJ VJ 1 . J 


0.0136 


(58.0) 




(0" 


0") 


6 /; (cm) 


10 p ' m 


13362+4831 fNGC 52561 


1 3 38 1 7 90 

1 J J O 1 / . 7V7 


4-48 16 41 

| lu 1 VJ 1 .VJ 


0.0278 


(120.0) 




(0" 


0") 


S'Vrm cm 1 

O VL.111, Alliy 


5 P 13 m ,6 m 


1 1428+5608 fMrk 2731 

1 _J ^ O l ^ V7 V J O 11V11JV ^ / ^7 / 


1 1 44 42 12 

1J I I Z. . 1Z 


4-55 53 1 3 5 

XJJ JJ 1 J.J 


0.0378 


(164.4) 




(0" 


0") 


3"(co) 


jp,m 


H470+3530* (\1GC 87391 

1JT/ WTJJJU ^ UVJV, o i -jy ) 


1 1 49 1 4 20 

ij ^xy ± ^ . z. \7 


4-35 1 5 21 

T-'J U Z.J.VJ 


0.0168 


(71.9) 


( 


-2.5" 


4-1 5") 


1 1 "(cm xl 

J_ _L V Vvlll,A / 


21 p 13 m 


F13500+3141 (3P7931 


11 52 17 87 
lj 1 / . 


-1-11 76 46 4 

TJl ^.VJ H-VJ.H- 


0.0446 


(194.9) 




(0" 


0"1 


i v>o; 


28 p ' m 


F1 3564+3741* (Nnr 53941 


1 1 58 11 60 

ID JO JJ.\J\) 


-U17 77 1 1 
j / z, / i j.vj 


0.0125 


(53.3) 




(0" 


0"1 
u ) 


5 /; (cm) 


21 p 13 m 


14003+3745* l"Nfir 54331 


1 4 07 16 00 

1+ UZ. JU.VAJ 


-1-17 10 18 


0.0145 


(61.9) 




(0" 


0"1 
u ) 


O.O (_Clliy 


21 p 13 m 


14151+2705* (Mrk6731 

y^y^jy \ t vj^j \ ivii jv vj / *j / 


14 17 21 00 

1^ 1 / Z. 1 . V7V7 


j_76 S 1 78 

^^Z.VJ J 1 Z.O.VJ 


0.0366 


(159.0) 




(0" 


0"1 




5P,m 


1 4 1 78-1-4Q77* ("7w747 070 Mrt 1 4Q0^ 


14 1 Q 41 70 
1^ iy H-J.Z.U 


-U4Q 14 17 


0.0256 


(110.3) 




(0" 


ft"\ 
u ) 


2" (cm) 


21 p 13 m 


14280+3126* (NGC5653) 


14 30 10.40 


+31 12 54.0 


0.0119 


(50.7) 


( 


-2.0" 


+1.6") 


17" (cm) 


21 p , 13 m 


14348-1447 


14 37 38.32 


-15 00 22.7 


0.0825 


(370.7) 




(0" 


0") 


< 2",3.5"(co,dbl) 


29 p . m I3 m 


15107+0724 (Zw049.057) 


15 13 13.07 


+07 13 32.0 


0.0129 


(55.0) 




(0" 


0") 


5"(co) 


23 p,m 


15163+4255* (Mrk848, Zwl07) 


15 18 06.20 


+42 44 42.0 


0.0402 


(175.1) 




(0" 


0") 


7"(sm,cm) 


5 p ,6 m , 13 m 


15243+4150* (NGC 5930, Arp 090) 


15 26 07.90 


+41 40 34.0 


0.0089 


(37.8) 




(0" 


0") 


2. 5" (cm) 


21 p ,13 m 


15322+1521* (NGC5953p, Arp091) 


15 34 32.30 


+ 15 11 38.0 


0.0065 


(27.6) 




(0" 


0") 


10"(cm) 


21 p ,13 m 


15327+2340 (Arp 220) 


15 34 57.24 


+23 30 11.2 


0.0182 


(78.0) 




(0" 


0") 


1.8"(co) 


2P,m 


15437+0234* (NGC 5990) 


15 46 16.50 


+02 24 56.0 


0.0128 


(54.6) 


( 


-2.1" 


-0.8") 


ll"(cm) 


21 p ,13 m 


16104+5235* (NGC 6090, Mrk496) 


16 11 40.70 


+52 27 25.0 


0.0292 


(126.1) 




(0" 


0") 


6" (cm) 


5 p ,13 m 


16284+0411* (MCG +01-42-008) 


16 30 56.50 


+04 04 59.0 


0.0245 


(105.5) 




(0" 


0") 


3. 9" (cm) 


21 p ,13 m 


16504+0228 (NGC 6240) 


16 52 59.05 


+02 24 05.8 


0.0243 


(104.6) 




(0" 


0") 


3" (co) 


7 p ,30 m 


17132+5313 


17 14 20.48 


+53 10 31.4 


0.0507 


(222.6) 




(0",- 


-1.0") 


2. 4" (cm) 


31 p ,13 m 


17208-0014 


17 23 21.92 


-00 17 00.7 


0.0428 


(186.8) 




(0" 


0") 


1.7"(co) 


2P,m 
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Table 3 - continued Observational parameters of the combined LIRG sample 



Name" 


RA (J2000) b 


Dec (J2000) 6 




(Ae a ,A8 s ) d 




Refs-f 


18425+6036* (NGC6701) 


18 43 12.27 


+60 39 10.5 


0.0132 (56.3) 


(+1.5", +2.0") 


15" (cm) 


9 p ,10 m 


19458+0944 


19 48 15.47 


+09 52 01.3 


0.1000(454.8) 


(0",0") 


< 0.8" (ir) 


14P,16 m 


20550+1656* (IIZw96) 


20 57 23.70 


+ 17 07 44.0 


0.0363 (157.7) 


(0",0") 


20"(sm,x) 


5P,6 m 


22491-1808 


22 51 49.86 


-17 52 24.4 


0.0773 (346.0) 


(-7.3", 0") 


2.5"(ir,cm,x) 


9 p ,15 m ,13 m 


23007+0836 (NGC7469) 


23 03 15.60 


+08 52 26.3 


0.0163 (69.7) 


(0",0") 


8" (co) 


32 p > m 


23365+3604 


23 39 01.25 


+36 21 08.4 


0.0644 (285.6) 


(0",0") 


0.95"(co) 


jp,m 



a IRAS name and the most common alternative(s), the asterisk marks sources for which CO line fluxes were obtained from an extensive 
literature search and line flux rectification process (see sections 2, 4.2). 
^Source coordinates (Rbeam) use d for CO observations. 

c The redshift used for receiver tuning (or the z co reported in the literature for LIRGs not in the original sample), and the corresponding 
luminosity distance in Mpc. 

d Offsets between expected CO source center R C m (assumed coincident with peak cm continuum) and observed position: 
AR=R cm — Rbeam (AR=0 when available CO or sub-mm images defined the CO source center, see 4.2). Thus R cm 
marks the true CO source center whenever AR^O. 

e CO region angular size (0s)=(0minfmaj) ' 5 , from interferometric maps. If these were not available then cm, near-IR, or sub-mm 
continuum images were used to set upper limits on (9 B ). The qualifiers for the images used are: (co), (cm), (ir), (sm) for CO, cm, IR, 
and sub-mm images, (x)=complex source morphology ((0 B ) then denotes the overall source size), (dbl)=double CO-bright nuclei. In 
the latter case (8 B ) refers to the largest one, and the second number denotes their separation. 

References used for source position (=p) and morphology/size (=m) information: l=Downes & Solomon 1998; 2=Zhu et al. 2009; 
3=Schinnerer et al. 1998; 4=Eckart et al. 1994; 5=Leech et al. 2010; 6=Mortier et al. 2009; 7=Wang et al. 1991; 8=Casoli et al. 
1992; 9=Sanders et al. 1991; 10=Condon et al. 1996; ll=Appleton et al. 2002; 12=Papadopoulos & Seaquist 1998a; 13=Condon 
et al. 1990; 14=Solomon et al. 1997; 15=Scoville et al. 2000; 16=Murphy et al. 1996; 17=Bryant & Scoville 1999; 18=Evans et 
al. 2002; 19=Sanders et al. 1988; 20=Wilson et al. 2008; 21=Yao et al. 2003; 22=Zenner & Lenzen 1993; 23=Planesas et al. 1991; 
24=Gracia-Carpio et al. 2007; 25=Evans et al. 2001; 26=Lisenfeld et al. 2000; 27=Aalto et al. 1997; 28=Evans et al. 1999; 29=Evans et 
al. 2000; 30=Tacconi et al. 1999; 31=Young et al. 1995; 32=Davies et al. 2004;33=Crawford et al. 1996 
9 Double nuclei in near-IR (Scoville et al. 2000), but only one is CO-bright (Evans et al. 2002). 
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Figure 4. The CO J=l-0 spectrum of IRAS 05189-2524, centered at 
cz=128()0kms~ 1 , obtained with the 12m telescope (see section 4.2.1). 



Indicatively, even for the ULIRG/QSO Mrk231 where ground and 
Herschel/FTS observations completed its CO SLED from J=l-0 up 
to J=13— 12 (van der Werf et al. 2010) it is always rj + u<rj j_i. 



5.1 Very high excitation gas: "out-exciting" Orion A and B 

The qualitatively different state of the molecular gas in (U)LIRGs 
than in spiral disks, has been noted in the past by several groups 
(e.g. Aalto et al. 1995; Solomon et al. 1997; Downes & Solomon 
1998; Papadopoulos & Seaquist 1998; Yao et al. 2003). Moreover 
ISM conditions similar to those in ULIRGs have been revealed in 
the nuclei of nearby spirals but involving only small fractions of 
their total molecular gas (e.g. Bradford et al. 2003). Nevertheless 
most such studies lacked J=3-2 and higher-J CO lines, and did not 
benchmak their very limited CO SLEDs to any regions where exci- 
tation mechanisms are well-determined (e.g. SF hot-spots in Galac- 
tic GMCs). 

For the highly excited CO SLEDs of some ULIRGs a super- 
position of two optically thick CO SLEDs, a high-excitation and a 
low-excitation one, can at most yield a "flattening" or a slower de- 
cline of high-J CO line ratios with respect to those expected from 
a low-excitation phase. Only dense and warm gas with at least par- 
tially optically thin CO lines can produce r43>r 3 2 and/or r32>r2i 
and the high excitation tails seen in Figure 5 (r2i,32>l)- Such CO 
line ratios are rare, and reported in the past only for the starburst 
nuclei of M 82 and NGC 3310 (Knapp et al. 1980; Olofsson & Ry- 
dbeck 1984; Braine et al. 1993), and although confirmed at lower 
values, they are still >1 (Weiss et al. 2001; Zhu et al. 2009), imply- 
ing extraordinary ISM conditions. Our study recovers this as a sta- 
tistically robust result for our sample. Indeed, while measurement 
errors still allow rj+i ; j<l for individual objects with measured 
r,j+i,j > 1 (J=l,2), this cannot be so for all the LIRGs in the high ex- 
citation tails in Figure 5 as the probability for this is: P2i~5x 10~ 7 
and P 32 ~4.4x 10~ 7 for the (2-l)/(l-0) and (3-2)/(l-0) line ratios 
(assuming independent gaussian probability distributions for each 
measurement). Thus some extraordinary physical conditions are 
possible for the bulk of the molecular gas of some (U)LIRGs. 
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Table 4. 12 C0 J=l-0, 2-1, 3-2 data for the combined LIRG sample 



Name 


CO J=l-O a 


CO J=2-l a 


CO J=3-2 a 


Refs. b 


f\f\r\ci ■ a no 1 

0005 / +402 1 


46 ± 7 


178 ± M 


2o7 ± 45(1.08/, Cj) 


x,l 


nnioo AO/in* 

UUJ22— 0840 


rr a a I 1 <it:c 

500 ± 125 


AAAA 1 CTAAC 

2000 ± 500 


of i a I r* cr ac 

2610 ± o50 


2 


nncnn i 1 OOC 
00509+1225 


o4 ± /(1. 11, K c ) 


"11/1 _L OQ/1 11 IT \ 

114 ± 2o(l.ll, K c ) 


Qr/; j_ 1 no/1 o"7 T/" A 

o50 ± 10o(1.2/, K c ) 


x,3 


m nci 1 T/i a * 
01053—1 /40 


A "1 _L 1 Q O 

byi ± loo 




O QO A _L aa A 
oo24 ± 000 


A C 

4,5 


n i (Yin 1 7fn* 
010/ /—I /O/ 


1 ( 5 it oO 




C^AV _L 1 AA 
50 ( ± 100 


A C 

4,5 


m /i 1 o i i Ac 1 * 
01418+105 1 


75 ± 15 




417 ± 83 


A C 

4,5 


020/ 1+385 / 


A AO _l_ Q Q 

40o ± OO 


"1 O Q Q _L O/lod 

12o5 ± 245 


2500 ± 515 


x,0 


CIO AO A i "3 TIC* 

02080+3 725 


147 ± 29 




r O/^ 1 1 1 C / 1 A CO T^" \ 

5oo ± 115(1. 052, K s h) 


6,7 


AT 11/1 _i_fl/1 
021 14+0450 


1 Kfl J_ QO 

loy ± o2 




1 1 aa. i o a a 
liyO ± 240 


A C 
4,5 


aoqo 1 adao, 


/1A _|_ a 


ZlO it 00(1. 30/, rs^c) 


O/l Pi -U c;^/0 AST f ^ e 
Z40 it O0(Z.U0/, lv c ) 


v on 


m /i n 1 aa n* 
02401—0013 


noon _L /I O/l 
ZooO ± 424 


( 1 1 Q _L A QO\ \/ 1 a4 

(1.10 ±0.22) X 10 


(1.(1 it O.o4) X 10 


in/ 
10 J 


at /i o i i /i q ao 
024 o o +4 302 


30 ± 4 


116 ± 21 


1 1 A _L 1 T{\ fl07 

liy ± 1 /(1. 08 /, Ij) 


«- 1 A 1 1 

X, 1,4, 1 1 


AOC1 Oj_1 /I A A* 

02512+1440 


1 CIA _l_ QO 

iy4 ± oy 




QQ7 _|_ 77 

o5( it / / 


/I c 
4, 5 


03359+1523 


155 ± 2o y 




256 ± 51 


4,7 


A/i o i o _i_ 1 /iiA 
04232+1430 


Q/l J_ T 

o4 it / 




Qor _L ia 

oo5 it /0 


V 1 1 

X, 1 1 


05083+ /930 


87 ± 9 


252 ± oo 


O/iA J_ 1 /ic:/i AC7 
o40 ± 145(1.08/, Lj) 


^-1 11 n 
X,l,l 1,12 


nc i on 

05189—2524 


4o ± / 


1 on _L OA/1 1A T/" \ 

loO ± 22(1.10, K s h) 


nci; J_ 0£V1 17 1/" \ 

25b ± ob(1.2/, K s h) 


X 


AO AO. A i CO /I Q 

08030+5243 


30 ± 6 


77 I 1 O 

(7 ± 12 


AO _L 1 'Z/l AO"7 /~*\ 

yo ± 1 ((1.08 /, Cj) 


1 1 

x, 1 1 


AO 1C/1 _lO C C C * 
08354+2555 


1 ftO J_ OA 

102 ± 20 


OfiT 1 -1 AZ 

20 ( ± 40 


£11 /I J_ AO 

014 it y2 


"7 11 1/1 1 C 
/, 13,14,15 


fiocao , on 1 C 

085/2+3915 


10.5 ± 1.5 


,1 1 I 1 O 

41 ± 12 




„ 1 1 1 A 

x,l 1,10 


no 1 oa_i_/1 /in 
09 120+4432 


1 fiO _L OQ 

102 ± 20 


/ /y ± 120(3.08, lv c ) J 


77ir _L i on 
t ID ± loO 


v A 1 1 A 1 "7 
X,4, /, 14,1 / 


(iW nucleus) 






540 ± 61 


x, 4 


" " (NE nucleus) 






oo c; _L A A 
2o5 ± 40 


x, 4 


AA1 A i A 1 Q /I 

0y320+0134 


70 ± 14 


345 ± 62 


ron 1 AFC/I AO"7 

5©y ± y5(i.08/, u) 


x,l 1 


AflCOt i 1 AAA* 

0y580+1000 


1 Oa _L OG/1 AOA f" \ 

120 ± 20(1. 080, K s h) 




007 J_ 1 OC/1 AOA T/" \ 

oot ± 120(1. 080, Kgh) 


x,0,18 


1AAQA nio* 
1003y— 3338 


65 ± 14 




316 ± 65 


A C 

4, 5 


i nmc i /i o5i 
10035+4852 


A o 1 in/1 OA iv- \ 

4o ± 10(1.20, K- S h) 


-| _L OA/1 OA IT \ 

loo ± 2\){l.ZU, Kg^) 




„ 1 1 
x, 1 1 


1 A1 T2 i AOOO 

101 /3+0828 


56 ± 10 




107 _L O/l /1 AO"7 /^\ 

127 ± 24 (1.08/, Cj) 


x, /,18 


1 A1 QA_i_1 no 

ioiyo+i 322 


□ ( ± ( 


1 AA _L OA/1 AA Ij^ \ 
100 ± 20(1.00, l^c) 


not; _L Q7 
220 it oi 


X,4, 19 


1 AQCA_i_C1/1 C* 
10350+5345 


1 AC\ _l_ OQC 

140 ± 25 


QOO _L 1 fifi 

522 ± 100 


1 coc; J_ QAAC 

1520 it oOO 




2 


i nc/cc i o /i /i o 
10505+2448 


/ 7 ± o 


o 07 _L o c; 
o27 ± o5 


500 ± 50 


^,1 /1"711 IO OA 

X,l,4, /,! 1,18,20 


1 1 1 A 1 i 1 1 AA 

i i iyi+1200 


4.5 ± 0.8 




OO _L CT / 1 A/1C f \ 

22 ib 5(1.045, K c ) 


X,21 


1 1 Tl 1 _i_1 A CA* 

1 1231+1450 


1 /IF 1 OA 
14o ± 20 




/1 1 Q _L Q7/ 1 A"7 Ti' ~\ 
410 it 5/(1.0/, 1V C ) 


7 10 


1 1 i COCA 

1 125 /+5850 


c. Q£; _L 1 1 c: 
5o0 ± 115 




4o00 ± 055 


«- A C 1 C OO 

X, 4, 5, 15,22 


1 OAA1 _i_A0 1 C* 

1 200 1 +02 1 5 


loo ± 25 




oil it 05 


1 Q 

18 


nil OiAIAC 

121 12+0305 


/lO J_ K 

42 ± 5 


152 ± 00(1.45, IVgh) 


KftQ I 11 Q/1 OA Ti^ \ 

5y5 it 110(1.20, IVgh) 


v "7 1 A 
X, /,10 


1 m/i aao/1 * 
12224—0024 


OO _L A 

22 ± y 




OQtr j_ co 
205 it 52 


1 Q 
18 


100/11 AA1A* 

12243—0030 


1 Q O J_ OQ 

lo2 ± 25 




not; j_ OA/1 

yy5 ± 204 


1 
l 8 


1 OC/1 A i C7AO 

12540+5 /08 


88 ± 9 


315 ± 30 


CCO _L OA 

50o ± oO 


OQ 

23 


1 0. AA 1 IQQfl* 

13001— 233y 


152 ± 30 




Qa a _i_ 1 7*0 
o04 ± 1 ( 


/I c 

4,5 


1 1 i ao i nci * 

13102+1251 


94 ± 19 




A O A 1 AA 

424 ± 90 


18 


Arp 238* 


110 ± 22 




onn _L cn 

oOO ± 00 


A C 

4, 5 


10.101 i 1 /I OT. 

13183+3423 


194 ± 16 


OCA _L 1 OA 

550 ± loO 


1 OA A _l_ 1 7*1 / 1 AO "7 

12y4 ± 1/1(1. 08 /, Cj) 


_ 1 7 11 111QO/10C 

X,l, /,! 1,13, 18,24,25 


111 QQ_i_AAQA* 

13188+0030 


1 ca _l qa 
150 ± oO 




VAfi 1 1 10 

(00 it 140 


1 Q 

18 


i no A OA1 /I 

1322y— 2y34 


Q Q O _L A Q 

052 ± 45 


1 OOC _l_ 1 OA 

12o0 ± 120 


1 Ai^A _L OAcfc 

iyoo ± 2y5 


OA 

X,20 


no *;o i /I Ol 1 

13302+4831 


175 ± 20 




A1 _i_ c\c\k 

912 ± 90 


„ A C 1 A 1 "7 

x, 4, 5, 14,17 


1 0./1OQ i CAAO 

13428+5008 


82 ± 9 


070 1 oc 

270 ± o5 


/ion 1 00/1 no*7 r~"i 
4o2 ± 52(1.08/, Cj) 


« 1 "7 11 

X,l, /,! 1 


1 Q/ITA i 1 CIA* 

134/0+3530 


1 1 O _l_ 0/1/111 f" \ 

llo ± 24(1.1 1, K s h) 




110c _L OQ1/1 11 f ^ 

1120 ± 2ol(l.l l,K s h) 


1 
18 


rnc aa i 1 1 a 1 
r 13500+3141 


52 ± 8 


1 re I OC/1 A1C V \ 

155 ± 25(1.035, K. c ) 


OAO _l_ CC/1 AA f \ 

20o ± 55(1.0y,R c ) 


x,27 


d 1 C A/1 117/11* 

r 13504+3 /41 


244 ± 51 




7CC J_ 1 cn 

/50 ± 150 


1 
1 8 


14003+3245* 


58 ±17 




853 ± 181 


18 


14151+2705* 


78 ±16 




340 ± 70 


4 


14178+4927* 


68 ±15 




544 ± 115 


18 


14280+3126* 


187 ±27 




806 ± 162 


7,18 


14348-1447 


56 ±8 


212 ±35 


360 ± 67 


x,7,28,29 


15107+0724 


120 ± 11 


605 ± 92 


710 ± 100(1.045, K c ) 


x,7,14,18,30 


15163+4255* 


114 ±23 




364 ± 73 


4,5 


15243+4150* 


52 ±14 




246 ± 55 


18 


NGC5953p* 


110 ± 24 




813 ± 166 


18 


15327+2340 


419 ± 36 


1127 ±69 


3674 ± 405 


31 


15437+0234* 


224 ±52(1.077, K sh ) 




636 ± 133(1.077, K sh ) 


18 


16104+5235* 


155 ± 31 




1062 ± 212 


4,5 


16284+0411* 


63 ±15 




580 ± 120 


18 
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Table 4 - continued 12 CO J=l-0, 2-1, 3-2 data for the combined LIRG sample 



Name 


CO J=l-0 a 


CO J=2-l a 


CO J=3-2 a 


Refs. 6 


16504+0228 


322 ± 29 


1492 ± 253 


3205 ± 642 


31 


17132+5313 


127 ±18 


168 ± 29 


264 ± 48 


x,14,17 


17208-0014 


160 ± 16 


688 ± 109 


1198 ± 190(1.087, G) 


x, 1,1 1,29,30 


18425+6036* 


234 ± 47 




1440 ± 288 


6,7 


19458+0944 


29 ±6 




261 ±65(1.087, G) 


x,ll 


20550+1656* 


121 ± 24 




527 ± 105 


4,5 


22491-1808 


33 ±6 


145 ± 34(1.36, K sh ) 


558 ± 190(2.1, K sh ) ( 


x,7 


23007+0836 


298 ± 27 


890 ± 103 


1600 ±240(1.12, K c ) 


x,7,13,14,26,32,33,34 


23365+3604 


39 ±6 


117 ±20 


288 ±71(1.087, G) 


x,l,ll 



"Velocity-integrated line flux densities in Jykms , with the value and type of any applied cor- 
rections reported in the parentheses (see 4.1), G: the pointing error bias, K sll =K sll (A0) (with 

A0=[(A(? Q ) 2 + (Adg ) 2 ] ): position offset correction, K c : beam-source geometric coupling cor- 
rection for sources with known CO or sub-mm dust emission source size 9 B . 

ft x=this work, l=Downes & Solomon 1998; 2=Zhu et al. 2009; 3=Barvainis et al. 1989; 4=Leech et al. 2010; 
5=Gao et al. 1999; 6=Narayanan et al. 2005; 7=Sanders et al. 1991; 8=Alloin et al. 1992; 9=Appleton et al. 
2002; 10=Papadopoulos & Seaquist 1998a; ll=Solomon et al. 1997; 12=Scoville et al. 1989; 13=Bryant & 
Scoville 1999; 14= Young et al. 1995; 15=Wilson et al. 2008; 16=Evans et al. 2002; 17=Sanders et al. 1986; 
18=Yao et al. 2003; 19=Gracia-Carpio et al. 2007; 20=Bayet et al. 2006; 21=Evans et al. 2001; 22=Casoli et 
al. 1992; 23=Papadopoulos et al. 2007; 24=Evans et al. 2005; 25=Mazzarela et al. 1993; 26=Papadopoulos 
& Seaquist 1998b; 27=Evans et al. 1999; 28=Evans et al. 2000; 29=Mirabel et al. 1990; 30=Planesas et al. 
1991; 31=Greve et al. 2009; 32=Davies et al. 2004; 33=Meixner et al. 1990; 34=Papadopoulos & Allen 2002 
c J=l-0, 3-2 for the central 40", J=2-l for the inner 20" (from Zhu et al. 2009). 
d From CO J=l-0 and the ((2-l)/(l-0)> ratio over the inner 22" (Casoli et al. 1992). 
e Probably an underestimate (source significantly larger than the JCMT beam). 
f Fluxes for the inner 20" of the galaxy. 

9 Shift factor K sh =1.25 applied for the Sanders et al. 1991 value. 
h Za upper limit. 

! Flux from an SMA image most likely collecting all of it given its single dish CO 3-2 flux is also recovered 
by the interferometer. 

J Uncertain value, large correction of a single CO 2-1 observation of a widely separated pair, pointing ~6" 
away from SW nucleus. 

k A 3x3 grid map at 7" grid size was made to obtain total line flux. 
; A large correction factor makes this value uncertain. 



In the Galaxy such high CO line ratios are found only in iso- 
lated "hot" spots in GMCs as results of strong and very localized 
excitation mechanisms. These are intense far-UV radiation from 
0,B stars inducing Photon-dominated Regions (PDRs) contain- 
ing warm gas near HII regions (Sakamoto et al. 1994), or super- 
novae remant (SNR)-induced shocks (e.g. Seta et al. 1998; Arikawa 
et al. 1999; Bolatto et al. 2003) which yield warm and strongly 
kinematically stirred gas (K v ; r ^>l). In the Orion A, B molecular 
clouds, the site of active formation of O, B stars in the Galaxy, the 
cloud-average (r2i)A,B~0.62-0.75, with r2i~0.90-1.3 found only 
in few isolated PDRs near HII regions and O, B star associations 
(Sakamoto et al. 1994). Thus r 2 i>0.8, observed in several LIRGs, 
already places their entire molecular gas reservoirs higher than the 
Orion A, and B clouds in terms of average excitation, while r2i > 1 
puts them on par with their "hottest" SF spots and the excitation 
state of molecular gas in SNR-GMC interfaces! This is indeed ex- 
traordinary given that the extreme conditions needed to produce op- 
tically thin CO lines in LIRGs must then be maintained over large 
molecular gas reservoirs (~(10 9 -10 10 ) M Q ). The CO/ 13 CO J=l-0, 
2-1 line ratios (Figure 7), with (Rio}~(R,2i}=18, are significantly 
larger than typical GMCs in the Galactic disk (Rio~4-6, e.g. 
Solomon et al. 1979; Polk et al. 1988), and independently indicate 
low optical depths. Using again the Orion A and B clouds to indi- 
cate how extraordinary such large Rio, 21 ratios are for galaxy-sized 
molecular gas reservoirs we note that R2i=10 is observed only in 



"hot" isolated interfaces of molecular gas to 0,B stellar associa- 
tions while for the bulk of the gas in Orion R2i~4-6 (Sakamoto et 
al. 1994). Even in the highly turbulent and warm molecular clouds 
in the Galactic Center the CO line emission remains mostly op- 
tically thick with Rio~5 and n ~5 (Oka et al. 1998). In other 
galactic nuclei such a dense and warm gas phase amounts to only 
~(l-2)% of the total molecular gas mass of a typical spiral disk. 

The extreme thermal and/or kinematic conditions needed to 
achieve even partially optically thin CO SLEDs can be simply 
demonstrated using the CO J=l— optical depth: 

-E 10 /k B T k - 



Tlu 



giA 1( 



1 - e" 



CO 



H 2 



n(H 2 ) 
(dV/dr 



(8) 



(assuming LTE), where Zlte is the partition function. After sub- 
stituting the various physical constants this becomes 



no = 1.29 x 10" 



-5.5/T1 



T k 



CO 



H 2 



n(H 2 ) 
(dV/dr) ' 



(9) 
4 (for 



where dV/dr is in units of kms -1 pc _1 . For [CO/H2]~10 
Solar metallicities), and setting the average velocity gradient of the 
gas phase as (dV/dr)=K v i r x (dV/dr) v i r where 



(-) 

V dr / vir 



n(H 2 



10 3 cm- 3 



1/2 



kms 1 pc 1 , 



(10) 



is the average velocity gradient expected for a virialized gas phase 
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Table 5. The CO J=4-3, 6-5 data 



Name 


CO J=4-3 a 


CO J=6-5 a 


00057+4021 


1049 ± 370 (1.15, G) 


< 380 6 (1 29 Gl 


00509+1225 




430 ± 155 (1.97, G) c 


02483+4302 




401 ± 131 (1.27, G) 


04232+1436 




< 3456 (1.21,G) 


05083+7936 


2033 ± 593 (1.15, G) 


688 ± 207 (1.22, G) 


05189-2524 


< 689 6 (1 47 K v.) 


638 ± 220(1.21,G) 


08572+3915 (NW nucleus) 




405 ± 158(1.21,G) 


09126+4432 CNF nucleus"! 




705 + 230 fl 90 K u) d 


09320+6134 




< 422 6 fl 35 Gl 


10173+0828 




< 156 6 (1 27 Gl 


10565+2448 




506 ± 148 (1.35, G) 


1 1257+5850 


6340 + 1 5S5 e 




11191+1200 




261 ±82 (1.14, K c ) 


12112+0305 (NE nucleus) 




327 ±143^ (1.22, G) 


12112+0305 (SW nucleus) 




< 123 6 (1.21.G) 


12540+5708 


1127 ± 265 (1.15, G) 


1320 ±400(1.29, G) 


13183+3423 




< 630 6 (1.28, G) 


13428+5608 




488 ± 156(1.45,G) 


F13500+3141 


842 ±245 (1.15, K c ) 


1086 ±347 (1.29, K c ) 


15107+0724 


1122 ±262 (1.073, K c ) 


904 ±241 (1.14, K c ) 


15237+2340 




3130 ±810 (1.21, G) 


16504+0228 




3321 ± 860 (1.21.G) 


17208-0014 


2312 ±654 (1.15, G) 


340 ± 140 s (1.17, G) 


23007+0836 


3970 ±820 (1.19, K c ) 


2355 ± 590 


23365+3604 


< 489 b (1.15, G) 


295 ± 110(1.37,G)» 



"Velocity-integrated line flux densities in Jykms -1 , with the value and type of any 
corrections applied (section 4.1) reported in the parentheses. A few CO J=6-5 fluxes 
differ slightly from those in PI 0b, a result of additional data. 
b 3cr upper limit. 

c The high velocity part of a double horn line tentatively detected, large correction 
applied because of extended source, uncertain total line flux. 

d Possible detection, highly uncertain because a pointing offset of ~3.85" results to 
a beam-shift correction K s j 1 =1.90. 

e From a partially completed CO J=4-3 map, and thus a likely underestimate. 
f Tentative detection of the NE nucleus of a double source. 
9 Very uncertain value. 



(with a~0.5-2.5 depending on the assumed cloud density profile, 
Greve et al. 2009; Bryant & Scoville 1996), we obtain 



tio = 63a 



-1/2 



-5.5/T k 



T k 



K vir 



(11) 



For most gas in typical GMCs: T kin ~(15-30) K, n(H 2 )~(500- 
10 3 )cm~ 3 , and K vir ~l (virial gas motions), yielding ti () ~6.5-33 
(for a=l.5). Only high temperatures (T k i n >100K) and/or highly 
non-virial gas motions (K v i r >10) can produce tioS;1- Moreover 
the conditions necessary for rj+i,.j>l become even more ex- 
treme for higher-J lines since, along with the high temperatures 
needed (so that Ej+i,j/(kBT k )<Cl), such ratios need densities 
high enough to thermalize the J+l >J line while retaining moder- 
ate/low optical depths. The latter becomes progressively more dif- 
ficult for higher-J lines since rj+i i j~(J + 1) 2 tio (for LTE and 
Ej+i,j/(k B T k )<l). Thus to keep ti <1 (so that e.g. T43~16ti 
will remain <1 yielding an optically thin J=4-3 line) both high 
Tkin and highly non-virial gas motions are necessary. Indeed 
for T k =100K, Equation 11 yields t 10 ~2.8x 10 _2 K^ v /n(H 2 ), 
which for n(H 2 )>(l-2)x 10 4 cm -3 (the critical densities of the 
J=3-2, 4-3 lines) gives tio>(3-4)K^, which can be brought into 



the optically thin regime provided that K v i r >10, which corre- 
sponds to highly unbound dynamical states for the gas. 

As the star-forming Orion A and B molecular clouds demon- 
strate, with low CO line excitation and large optical depths prevail- 
ing when averaged over their entire volume, such extreme condi- 
tions cannot easily dominate the CO SLEDs of individual SF clouds 
let alone the global CO SLEDs of entire galaxies. Thus rj+i,j>l 
(J= 1,2,3) ratios are not expected, even in vigorously star-forming 
LIRGs, and are indicative of strong galaxy-wide effects maintain- 
ing high gas temperatures as well as highly non-virial velocity 
fields up to densities n>10 4 cm~ 3 throughout their large molec- 
ular gas reservoirs. As we will argue in more detail later such a 
high-excitation gas phase is irreducible to an ensemble-average of 
individual SF molecular clouds. For the rest of LIRGs more ordi- 
nary excitation states of SF (r2i,32<0.9) and non-SF (r2i~0.5-0.7, 
r32~0.3) gas adequately account for their observed CO SLEDs. 

High resolution submm imaging with ALMA will reveal ex- 
treme CO line excitation and thus extraordinary ISM conditions 
much more frequently as the high CO ratios of SF-powered re- 
gions or AGN-excited gas will no longer be "watered-down" by 
lower-excitation cloud ensembles. Line ratio imaging can then re- 
veal deeply dust-enshrouded AGN in galactic centers where nei- 
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Figure 5. The distributions of the CO (2-l)/(l-0) (left) and (3-2)/(l-0) (right) global brightness temperature ratios for the LIRGs in our sample. The shaded 
area marks the line ratio distributions for sources with CO emission region angular sizes of ^15" (the HPBW of the JCMT at 345 GHz). The measured range 
of the (2— 1)/(1— 0) line ratio in the hot spots near HII regions or PDRs around O, B stars in the Orion A, and B molecular clouds (Sakamoto et al. 1994), and 
the expected range of the CO (3-2)/(l-0) ratio in these high-excitation regions are indicated by the horizontal bar. 



ther optical nor even IR lines are capable of doing so. An early 
such demonstration is the imaging of the central ~100pc of grand 
design spiral galaxy M51 which finds high CO line excitation 
with r 32 ~1.9, and possibly caused by an AGN-powered jet (Mat- 
sushita et al. 2004). 

5.2 High-excitation gas in LIRGs: "bottom"-stirred by SF 
feedback? 

A warm (Tki n >100 K) unbound (K v i r »l) phase concomitant with 
a less excited cooler and denser one is a feature of standard two- 
phase models in (U)LIRGs (e.g. Aalto et al. 1995; Papadopoulos 
& Seaquist 1999), and could in principle explain the very high- 
excitation CO SLEDs in our sample. In these systems a frequent 
non-convergence of global CO line ratios (mostly CO, 13 CO 1-0, 
2-1) to an average ISM state is attributed to an envelope-like diffuse 
(n~(10 2 -10 3 ) cm~ 3 ) non self-gravitating warm gas with high tur- 
bulent linewidths (and thus small/moderate rio), surrounding much 
denser self-gravitating gas where CO and even CO lines have 
large optical depths (Aalto et al. 1995). However, the difficulty of 
maintaining low optical depths and well-excited CO J=3-2, 4-3 
lines (so that 143,32 >1) makes this two-phase ISM model problem- 
atic as the "envelope" phase should then have n(H2)>10 4 cm~ 3 
(n crit (3-2,4-3)~(l-2)x 10 4 cm~ 3 ) rather than n<10 3 cm" 3 . 

SF-feedback effects capable of producing highly-excited CO 
SLEDs with moderate/low optical depths up to 1=3-2, 4-3 do occur 
in dense gas regions deep inside molecular clouds, and are driven 
by the massive stars. The warm PDRs and SNR-shocked gas re- 
gions, along with strong radiation pressure onto the dust grains 
mixed with the molecular gas can produce a warm and non self- 
gravitating (K v i r >l) dense phase near SF sites. It is worth noting 
that radiation pressure from massive stars is powerful enough to 
even be the main SF-regulator in the starbursts found in ULIRGs 
(Andrews & Thompson 201 1). Such a "bottom"-stirred ISM would 



yield a very different two-phase differentiation, with a warm, dense, 
and non-virial phase lying deeply embedded in molecular clouds, 
surrounded by a less dense cooler one extending further from the 
vicinity of SF sites. A warm, dense, and non-virial phase is indeed 
compatible with the CO line ratios of the SF "hot-spots" of Orion 
A,B molecular clouds (see 7.2.1), with CO SLEDs expected to be 
highly excited and partially optically thin up to high-J lines, a result 
of the high temperatures and K v i r values. 



5.2.1 A "bottom" -stirred ISM, some predictions 

In the context of the old 2-phase model, the high CO/ 13 CO 1=1-0 
ratios in ULIRGs (Rio>20) are attributed to the much more tur- 
bulent (K v i r ;^l) diffuse/warm "envelopes" than those of molec- 
ular clouds in less extreme star-forming systems (Rio~10-15) or 
SF-quiescent ones like the Milky Way (Rio~4-6) (Aalto et al. 
1995). A much larger amount of turbulent energy per gas mass 
is available in ULIRGs, a result of strong mergers, and its dissi- 
pation can readily warm molecular gas (Pan & Padoan 2009) and 
increase Rio. In a "bottom"-stirred ISM however, the dense phase 
is the one most affected by the SF feedback and thus high (>15) 
12 CO/ 13 CO Rj+i,j J+l=3, 4, 5, 6 andHCN/H 1:i CNJ=l-0 ra- 
tios are expected, as dense gas can now attain low/moderate molec- 
ular line optical depths. These transitions (which would be very 
optically thick for dense gas that has K v i r ~l), with n cr it~(10 4 - 
10 5 ) cm~ 3 , probe the dense gas with negligible contributions from 
diffuse cloud envelopes. High C/ 13 C-based isotopologue line ra- 
tios at ever increasing critical densities (e.g. HCN/H 13 CN J=3-2, 
n cr it~4 x 10 6 cm~ 3 ) would then provide an ever sharpening dis- 
tinction between a "bottom"-stirred and a standard two-phase ISM 
model in (U)LIRGs. 
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Table 6. The 13 


CO line data 






Name 


13 CO J=l-0 a 


13 CO J=2-l a 


References 6 


0005 /+4021 




<" O /I c 


• 

this work 


00309+1225 


4.0 ± 0.8 


(K=o ±2)" 


1, this work 


0105 J— 1 /4o 


/T) ~> -i v\c . d 
( 1 f ) ' 




L 


020 / / 


37 i 6 


(K=12 ± o) 


this work 


024U1— 00 i J 


loO ± 2o 


10o4 ± 20U 


J 


U24o3+4JU2 


<" o 1C 


<" QIC 

^ 21 


this work 


UoJjy+1 52:5 




/D 1 o _L Q\e 

(K=12 ± o) 


this work 


(jjuoJ+/yjo 


7.9 ± 2.3 


23.0 ± 4.5 


this work 


nc 1 on Tci/i 
05189-2524 




21 ± / 


this work 


08J54+2555 




/D 1 Q\C,d 


O 


09126+4432 




< 55 c 


this work 


nonet ■ <ci i >i 
09J20+O1 J4 




18 i 3 


this work 


1UU3j+4oj2 




o Q ^ o c; f 

y.y ± z.o j 


this work 


1017J+0828 




5.2 ± 1.0 


this work 


105o5+244o 


4.7 ± 1.3 


(K^ 18) 


this work, 6 


1 1257+5850 


(K=27 ± 5) 




7 


12540+5 /Ob 




6.1 ± 2.1 


this work 


13183+3423 




26 ±4 


this work 


13229-2934 


13.5 ±2.0 


87 ± 13 


this work, 8 


13428+5608 




36 ±9 


this work 


F13500+3141 


< 11. 8 C 




this work 


14348-1447 




16.7 ±5.6 


this work 


15107+0724 


6.8 ± 1.8 


23 ± 7 


this work 


15327+2340 


9±2 


59 ±7 


this work, 9 


16504+0228 


6.5 ± 1.9 


26 ±4 


this work, 9 


17208-0014 




< 18 c 


this work 


23007+0836 




54 ± 13 


8 


23365+3604 




6.2 ± 1.4 


this work 



"The velocity-integrated line flux densities in Jy km s 1 . 
b l=Eckart et al. 1994; 2=Aalto et al. 1995; 3=Aalto et al. 1991; 
4=Casoli et al. 1992; 5=Papadopoulos & Seaquist 1998a; 6=Wil- 
son et al. 2008; 7=Aalto et al. 1997; 8=Papadopoulos & Seaquist 
1998b; 9=Greve et al. 2009. 
c 3cr upper limit. 

d Only R= 12 CO/ 13 CO is available (reported in parenthesis) for 
areas smaller than the total CO-emitting region (see 4.2). 
e Large pointing offset for both CO and 13 CO line measurements 
makes absolute fluxes uncertain but leaves line ratio invariant. 
^Uncertain value, large K sh correction applied (K sh =2). 



5.3 The quest for a new ISM excitation mechanism in 
ULIRGs 

The small amount of gas mass typically found in the SF "hot- 
spots" of ordinary star-forming GMCs renders even a "bottom"- 
stirred ISM model problematic in producing the highest excitation 
CO SLEDs found in our sample. Indeed, even if all the dense gas 
(>10 4 cm- 3 )in GMCs was feedback-affected it would still amount 
only to ~few% of their total mass. This cannot "skew" their CO 
SLED anywhere close to the high excitation levels observed for 
some ULIRGs (e.g. IRAS 12112+0305, IRAS 17208-0014). Much 
larger mass fractions of GMCs must be dense and SF-active for 
such SLEDs to emerge. High cloud boundary pressures and tidal 
striping of outer GMC envelopes in a merger environment can pro- 
duce much denser clouds, a situation exemplified by the Galactic 
Center (Giisten & Phillip 1994). Large amounts of molecular gas 
(~10 10 M ) at densities ~(10 5 -10 6 ) cm~ 3 , and a possible steep- 
ening of the (average density)-size power law (n)ocL~ k to k>l 
(k=l for Galactic GMCs, Larson 1981) certainly argue for a much 
denser molecular ISM in ULIRGs (Greve et al. 2009). At the same 



time this hinders far-UV photon propagation even more (i.e. the 
PDRs will be smaller) while SNR-induced shocks in such dense 
ISM will dissipate fast and produce shock-excited molecular gas 
regions that are even more confined around the SNR sites. Thus a 
quest for mechanisms that can globally heat the large amounts of 
dense gas found in ULIRGs to high temperatures opens up. 



5.4 Strong AGN feedback: a new global driver of ISM 
excitation 

Two intriguing objects, the optically/IR-luminous QSO 
PG 11 19+120, and the powerful FRII radio galaxy 3C293, 
widen the range of possibilities regarding globally operating 
molecular gas excitation mechanisms in galaxies. This is partic- 
ularly clear in 3C 293 where a strong jet-ISM interaction injects 
large amounts of mechanical energy into its highly turbulent 
ISM and powers luminous CO J=4-3, 6-5 lines in an otherwise 
SF-quiescent ISM with subthermal low-I CO line ratios and 
cold dust emission (Papadopoulos et al. 2008, 2010b). Strong, 
penetrating irradiation by an X-ray luminous AGN can volumet- 
rically heat molecular gas by creating giant X-ray Disossiation 
Regions (XDRs), and is another powerful and possibly global 
ISM excitation mechanism in AGN-hosting LIRGs (Meijerink 
et al. 2007; Schleicher et al. 2010). Unfortunately there are no 
X-ray observations of PG 1 1 19+120 to assess the X-ray luminosity 
of its AGN, while we find no particular predominance of AGN 
amongst the LIRGs populating the high-excitation tails in Figure 
5. Very high-J CO SLEDs (I>8-10) will provide the tools for 
distinguishing between AGN-induced XDRs, SF-induced PDRs 
as the cause of high ISM excitation (Meijerink et al. 2006; van 
der Werf et al. 2010). Finally the recent discovery of a massive 
molecular gas outflow, likely driven by the AGN in Mrk231, by 
Aalto et al. (2012) inserts an exciting new element in this debate 
by demonstrating that much higher density gas (>10 4 cm~ 3 ) can 
also be affected by AGN feedback. 



6 CO LINE RADIATIVE TRANSFER MODELS 

Interpreting global CO line ratios in terms of underlying average 
gas physical conditions is the task of molecular line radiative trans- 
fer models. In the present analysis of the CO, 13 CO lines we will 
use a Large Velocity Gradient code (see Papadopoulos & Sequist 
1998 and references therein) that has been shown adept at discern- 
ing the average thermal and kinematic states of the molecular gas 
when constrained by CO, 13 CO line ratios (e.g. Mao et al. 2000; 
Weiss et al. 2001; Yao et al. 2003). Our LVG code uses the ob- 
served CO ratios R^ obs ' to find the regions in the [n(H2), Tk, A co ] 
parameter space where x 2 =Sk(l/fk) 2 [R^° bs ' ) — Rk] 2 (where Rk 
is a model line ratio and <Tk the measurement error) is minimum. 
A co =[CO/H 2 ]/(dV/dR) (in (kms^ 1 pc" 1 )" 1 ), where [CO/H 2 ] 
is the CO abundance and dV / dR is the average local velocity gra- 
dient of the turbulent CO line emitting gas. The temperature range 
we considered is: 15 K<Tki n < 150 K, with the minimum value typ- 
ical of cold SF-quiescent gas and the maximum of the Cold Neutral 
Medium HI. Moreover only LVG solutions with K v i r > 1 where 



(dV/dR) 
(dV/dR) 



1.54 



virial 



[CO/H 2 



n(H 2 ) 
10 3 cm- 3 



-1/2 



(12) 



(using equation 10) are accepted. The latter determines the aver- 
age dynamic state of the molecular gas, with K v i r ~l-3 for the 
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Table 7. 12 CO line ratios and SF-powered IR luminosities 



Name 


r 21 
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Table 7 - continued 12 CO line ratios and SF-powered IR luminosities 
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20550+1656* 




0.48 ±0.13 






7.10 ± 1.40 


4.71 (42) 


22491-1808 


1.10 ±0.32 


1.88±0.72' 1 






8.97 ± 1.63 


13.0 (49) 


23007+0836 


0.75 ±0.11 


0.60 ±0.10 


0.83 ±0.19 


0.22 ±0.06 


3.48 ±0.32 


1.61 (42) 


23365+3604 


0.75 ± 0.17 


0.82 ± 0.24 


< 0.78 d 


0.21 ± 0.09 h 


7.30 ± 1.12 


6.19 (61) 



"The CO line ratios r.I+l.I=Lco(.I+i-.I) / L CO(l-0) 

ft The CO 1-0 luminosity in units of LpKkms -1 pc 2 (see Equation 6). 

c The SF-powered part of L IH and the associated dust temperature (see PlOa for details). 

d 3-<x upper limit. 

e CO J=6-5 flux from the NW nucleus where most of the CO emission is (Evans et al. 2002). 
^Likely higher as CO 4-3 is obtained from a partially completed map. 
9 Estimated for the NE nucleus where 75% of CO(1-0) emission arises (Evans et al. 2002). 
^Uncertain value, see footnotes of Tables 4,5. 



Table 8. The CO/ 13 CO line ratios 



I 5 



Name 


J=l-() a 


J=2-l a 


00057+4021 




> 7 b 


00509+1225 


8±2 


6±2 


01053-1746* 


> yjb,c 




02071+3857 


10 ±2 


12 ±3 


02401-0013* 


14 ±3 


10 ±3 


02483+4302 


>^ b 


>5 b 


03359+1523 




12 ±3 


05083+7936 


10 ±3 


10 ±2 


05189-2524 




6±2 


08354+2555 




> l 8 b,c 


09126+4432 




>13 6 


09320+6134 




18 ±4 


10035+4852 




16 ±5 


10565+2448 


15 ±4 


> 18 6 


1 1257+5850 


27±5 C 




12540+5708 




47 ± 16 


13183+3423 




30 ±7 


13229-2934 


26 ± 5 


13 ±2 


13428+5608 




7±2 


F13500+3141 


> 4 b 




14348-1447 




12 ±4 


15107+0724 


16 ±4 


24 ±8 


15327+2340 


43 ± 10 


17 ±2 


16504+0228 


45 ± 14 


53 ± 11 


17208-0014 




> 35 6 


23007+0836 




15 ±4 


23365+3604 




17 ±5 



a h' (CO)/L' ( 13 CO) line ratios. 

6 Ratio that corresponds to 3-cr upper 

limit on 13 CO flux. 

c Ratio measured over a smaller area 

than that of total CO J=l-0 emitting 

region. 
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Figure 6. The distribution of the CO(6-5)/(l-0) brightness tempera- 
ture ratios for the LIRGs in our sample. The optically luminous QSO 
PG 11 19+120 has the highest excitation with rg 5 ~1.6. Some of the low 
ratios (^0.3) may have been surpressed by a strong dust continuum becom- 
ing optically thick even at short submm wavelengths (see PlOa). 



(mostly) self-gravitating GMCs in the Galactic disk (and K v i r <Cl 
corresponding to dynamically unattainable gas motions). There is 
no theoretical upper limit for K v i r , and K v i r ^>l in regions where 
strong injection of turbulent energy occurs (e.g. in GMC/SNR in- 
terfaces). Indicatively, for average densities of ~(10 2 -10 4 )cm -3 
K vir =20 corresponds to (dV/dR)~(5~50)kms _1 pc -1 (for a = 
1.5), encompassing the values observed in the highly turbu- 
lent gas of the Galactic Center (Dahmen et al. 1998; Oka et 
al. 1999), in pre-star-forming molecular clouds in the Galaxy 
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Figure 7. The joint distribution of the CO/ 13 CO J=l-0, 2-1 line ratios for 
the LIRGs in our sample. 



(Falgarone et al. 1998) and galactic SF nuclei where typically 
(dV/dR)~(2-6)kms _1 pc _1 (Weiss et al. 2001). Furthermore 
we adopt [CO/ 13 CO]=50 as appropriate for the environments in 
LIRGs (e.g. Henkel et al. 1993) and measured also for the Galactic 
disk in the Solar vicinity (Wilson & Penzias 1993), while we set 
[CO/H2]=10~ 4 (needed for computing K v i r from Equation 12). 
Finally we note that some CO J=6-5 line luminosities may be di- 
minished by significant dust optical depths even at short submm 
wavelengths (P10a,b), while in a few cases (e.g. Arp 193) they may 
be underestimated because of substantial pointing offsets. Thus we 
do not use them as constraints in our LVG modeling of the entire 
sample. Finally, we only made LVG models when N>3 lines were 
available, except in the special cases where rj+i,j>l that indicate 
quite unique high-excitation regions of the LVG parameter space. 

In Figure 8 we show the T kin , T r (l-0) (CO J=l-0 brightness 
temperature), n(H2), and K v i r distributions as derived from the 
LVG radiative transfer models, including the degenerate or poor 
solutions. A general n(H2)-Tki n degeneracy exists when r.j + i i j<l 
and 13 CO lines are unavailable as constraints, with such ratios 
compatible with dense/cold and with diffuse/warm phases (e.g. 
for IRAS 14348-144: 10 4 cm~ 3 /20 K and 300 cm~ 3 /150 K). These 
nevertheless typically have different K v i r , which can be used to dis- 
criminate between such disparate LVG solution ranges once 13 CO 
lines become available. The few states with n(H2)>10 4 cm~ 3 are 
associated with LIRGs where rj+i^j^l (the high excitation tails in 
Figure 5), they also have T k i n >100 K, and involve large molecular 
gas masses of ~(4-6) x 10 9 M Q . These extreme systems are studied 
in more detail in our final paper (Papadopoulos et al. 2012). 

The wide K v j r distribution is a result of the large range of 
Rio, 21 ratios (Figure 7), with average dynamical states having 
K v ir~0.7-2 (~virial values) to K v i r >20 (strongly unbound states). 
This diversity is expected given that the galaxies in our sample in- 
clude isolated disk-dominated systems (where GMCs are mostly 
self-gravitating) and strongly evolving mergers (where GMCs can 
be fully disrupted by tidal fields), and this has been noted from 
early studies of such systems (e.g. Aalto et al. 1995). The distri- 



bution of the average brightness temperature of the CO J=l-0 line 
T r (l-0) is rather narrow, with an average of ~16 K and o"(T r )~4K 
(for the distribution of the good, non-degenerate LVG solutions). 
The Xg-d=T kin /Tdu S t and Xth=T ex (3-2)/T kin distributions (Fig- 
ure 9) indicate how much "warmer" (or "colder") a CO SLED is 
with respect to the SED of the concomitant dust reservoir, and the 
thermalization level of the J=3-2 transition. 

For most LVG solutions T kin /T dust >1, which is expected 
since photoelectric, and even more so turbulent or CR heating of 
the ISM affect the dust much less. Only at very high densities 
(>10 5 cm~ 3 for Galactic ISM environments) thermal equilibrium 
between gas and dust is attained (T k i n ~Td us t)- Some large ther- 
mal decouplings (Xg-d>3) exist, with one of the highest found in 
the powerful radio galaxy 3C 293 where a strong jet-ISM interac- 
tion drives a "hot" CO SLED while its dust SED remains cold, 
reflecting Milky- Way levels of star formation (PlOb). Thus turbu- 
lent heating (in this case AGN-driven) can heat galaxy-sized molec- 
ular gas reservoirs to high temperatures, even in the absence of 
strong SF-powered radiation fields, an important issue discussed 
further in section 7.1. The xth values range from very subther- 
mal (~0. 1-0.2), to fully thermalized (~1) in the high-excitation 
ISM of some ULIRGs. Low xth values on the other hand are 
obtained mostly for a high temperature (T k i n ~(60-150) K), dif- 
fuse (n(H2)~(10 2 -10 3 ) cm -3 ) phase. This can be typical in turbu- 
lent GMC envelopes, and may even dominate the CO J= 1—0, 2- 
1 line emission in ULIRGs (DS98), yet it cannot contain much 
of their molecular gas mass. Indeed in such galaxies with SF- 
powered IR luminosities of LJ R ~10 L0 (see PI 0b) the dense 
gas phase associated with SF sites will have a minimum mass 

of M S F(H 2 )~L^ ) / e g,*~ 4 x !0 9 M © ( for an Eddington-limited SF 
efficiency of e g ,*~25OL /M ). This is already ~2xM(HI+H 2 ) of 
gas-rich spirals (e.g. Young & Scoville 1991), the typical progeni- 
tors of ULIRGs (Sanders & Ishida 2004), thus leaving little room 
for any significant additional molecular gas mass in a diffuse phase 
in such systems. A very dense dense star-forming gas phase con- 
taining most of the molecular gas mass is indeed strongly suggested 
by HCN J=l-0 observations (Gao & Solomon 2004) and conclu- 
sively shown for individual ULIRGs using heavy rotor (e.g. HCN, 
CS) and 13 CO lines (Papadopoulos et al. 2007; Greve et al. 2009). 



6.1 Low-excitation CO SLEDs and the average ISM states 

Only few "cold" CO SLEDs are found in our sample, indicated 
by the low excitation "tail" of CO line ratios with r2i ^0.6 and/or 
r 32^0.30 (Figure 5). This is expected for what at the end is an 
IR-selected (and thus SFR-selected) galaxy sample, with the limi- 
tations of such samples in selecting large numbers of (cold-ISM)- 
dominated systems discussed thoroughly by Dunne et al. 2000. 
Furthermore, unlike the high excitation CO SLEDs that uniquely 
correspond to warm, dense (and often unbound) gas states, sub- 
thermally excited "cold" low-J CO SLEDs cannot be uniquely at- 
tributed to a quiescent ISM state. Indeed they could mark a cold 
(~(10-20)K), self-gravitating (K v i r ~l) phase of moderate den- 
sities (~10 cm~ 3 ), typical of SF-quiescent GMCs, but also a 
warmer (>30K), low-density (<10 3 cm~ 3 ), and highly unbound 
gas phase (K v i r ^> 1) more typical of turbulent cloud "envelopes" or 
intercloud gas in ULIRGs (Aalto et al. 1995; Downes & Solomon 
1998). This renders low-excitation CO J+l^J, J+l<3 SLEDs 
highly degenerate to the average state of the ISM in galaxies. This 
is obvious in individual systems such as 3C 293 and Arp 220 where 
subthermal low-J CO lines give no hint about their massive and 
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Figure 8. The n(H2), K v i r , Tki n , and T r (l-0) distributions obtained from ( 
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t-phase LVG radiative transfer models of the CO lines of the LIRG sample, 
solutions obtained for each system (see section 6). 



highly-excited gas reservoirs revealed only via CO J+l — >-J, J+1^4 
or heavy rotor (e.g. HCN, CS) lines (Papadopoulos et al. 2008; 
Greve et al. 2009). 

This degeneracy casts doubts on recent claims for Galactic- 
type ISM in distant gas-rich star-forming disks (Dannerbauer et 
al. 2009; Daddi et al. 2010) as current observations do not ex- 
tend beyond CO J=3-2 where such degeneracies diminish (see 
section 7.2.1). Among the galaxies with low-excitation CO line 
ratios in our sample this was possible for only two LIRGs 
(IRAS 03359+1523 and IRAS 05 189-2524) where available 13 CO 
line observations (see Table 8) played a key role in determining 
Galaxy-type ISM states as the most likely ones. Observations of 
15 CO lines of distant systems like the BzK galaxies will soon be- 
come possible with ALMA and will be decisive in breaking the 
aforementioned degeneracies of low-J CO SLEDs with Rio,2i~5- 
10 corresponding to SF-quiescent Galactic-type ISM conditions 
while Rio,2i>15 indicating vigorous star-forming and/or merger 
ISM environments. 



7 CO SLEDS OF LIRGS: THEIR EXCITATION RANGE, 
AND POWER SOURCES 

The range of normalized Lj^ j=Lj+i, j /Lir CO SLEDs as deter- 
mined from our study is shown in Figure 10. From this it becomes 
obvious that: a) beyond J=3-2 very large variations become possi- 
ble, b) some global SLEDs are on par with those of the Orion SF 
"hot-spots". The highest excitation SLED is AGN-driven (3C293) 
followed by those for extreme starbursts (e.g. IRAS 17208-0014, 
IRAS 12112+0305). We must note that the large L^ j range, 
spanning ~3-5 orders of magnitude beyond J=4-3, reflects both 
the degeneracy of the one-phase radiative transfer models as well 
as the wide range of ISM excitation found in our sample. Recent 
Herschel/FTS observations of some these (U)LIRGs fully deter- 
mine this higher-J part of their CO SLEDs confirm this picture, and 
will be presented in a forthcoming paper. Finally from Figure 10 it 
is obvious that a positive K-correction for redshifted CO line emis- 
sion from distant galaxies due to a rising Lj 1 ^ with J-level (which 
makes possible the detection of very distant systems in high-J lines 
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by nearly countering luminosity-diminishing with distance) can be 
guaranteed only up to the J=3-2, 4-3 transitions. 

7.1 Cosmic rays and turbulence: the dominant molecular gas 
heating mechanisms in ULIRGs? 

What is to be made of those global CO SLEDs that are so highly 
excited that they are on par or even surpass those expected for 
the Orion SF "hot-spots"? Moreover, maintaining high temper- 
atures (Tki n >100K) for metal-rich gas with average densities 
n(H2)>10 5 cm -3 , as implied by some of these SLEDs, is dif- 
ficult since cooling Aoc[n(H2)] 2 (Lequeux 2005), and is thus 
~10 4 -10 6 times stronger than for average densities of molecular 
clouds (~10 2 -10 3 cm~ 3 ). The galaxy-sized gas reservoirs (~(2- 
5)x 10 9 Mq) involved make such highly-excited CO SLEDs even 
harder to maintain energetically since, for their high average densi- 
ties, far-UV radiation fields will be reduced by factors of ~ 10 4 over 
distances of <0. 1 pc. The corresponding SF-powered PDRs will 
thus be very localized, involving only small fractions of a CMC's 
mass, as demonstrated by the SF "hot-spots" in the Orion A, B 
clouds, while the same is expected for SNR-shocked regions (i.e. 
the "mechanical" counterpart of SF feedback). 

Cosmic rays (CRs) and/or turbulence, already implicated as 
powerful heating sources of the molecular gas in Galactic Center 
(Yusef-Zadeh et al. 2007), and powering the high-excitation CO 
SLED of the starburst nucleus of NGC253 (Bradford et al. 2003; 
Hailey-Dunsheath et al. 2008), can provide strong volumetric heat- 
ing on the mass scales necessary to power the high-excitation CO 
SLEDs of ULIRGs. For CRs this has been recently demonstrated, 
with the large CR energy densities of UcR~[(few)x 10 3 ]UcR,Gai 
expected in such galaxies penetrating much deeper into dense 
molecular regions than far-UV, optical and even IR photons and 
heating their mass up to ~100K (Papadopoulos 2010). The very 
turbulent gas disks found in such galaxies (DS98) (er z ~(40- 
140) km s -1 versus er z ~(5-6)kms~ 1 for ordinary face-on spirals) 
and the dissipation of their supersonic velocity fields can consider- 
ably warm the molecular gas throughout their volume. Following 
Pan & Padoan 2009, the average turbulent heating rate is 

Tturb = M m H<e(x,t)), (13) 

where /i=2.35 is the mean molecular weight for molecular gas, n 
is the gas density, and (e)=(e(x, t))=l/2(v / 3u v ) 3 /L is the average 
turbulence dissipation rate per unit mass (assuming an isotropic ve- 
locity field). Using the av-L (velocity dispersion)-(cloud size) rela- 
tion 



(fef(£)' 



(14) 



(Larson 1981), which is now well-established for the Galaxy with 
nearly constant o- =1.2kms _1 and power-law index a=l/2 (Heyer 
& Brunt 2004). The pressure-dependence of the normalization, 
with P o /kB=10 4 Kcm~ 3 , can be derived from the virial theo- 
rem applied for clouds with boundary pressure P e (Chieze 1987; 
Elmegreen 1989). This pressure-dependence is very important as 
it predicts a higher normalization of the cr v -L relation for high- 
pressure environments. This has actually been observed in the 
Galactic Center (Giisten & Phillip 2004) and astonishingly also in 
a molecular gas disk in a distant ULIRG at z~2.3 (Swinbank et al. 
201 1). From Equations 13 and 14 after replacing the expression for 
(e), setting a=l/2, and substituting units we obtain 



/P e \ 3/4 

Tturb = 3.3 x 10~ 27 ncr o ,n \jT) L pc er g 



where a ,n=o"o/(kms _1 ) and L pc =L/(pc). On the other hand the 
(mainly) molecular line cooling for self-gravitating gas is well- 
approximated by the expression: 



Aii nc = 1.9 x 10~ 27 n 



erg cm s 



(16) 



(Papadopoulos 2010 and references therein), with most of it due 
to the CO SLEDs though other molecules (e.g. 13 CO) significantly 
contribute when n>10 4 cm -3 (Goldsmith 2001). Thus for thermal 
equillibrium: r tU rb=Aii nc , we obtain 



Tki, 



12a 



1/3 



Po/k E 



10 4 Kcm- 3 



1/4 



1/6 



K. 



(17) 



-3 -1 

cm s 



(15) 



Hence, for the typical pressures expected in the Galaxy and the 
midplane of ordinary spirals (~10 4 Kcm~ 3 ), turbulent heating can 
easily account for the temperatures of UV- shielded regions (~(10- 
15) K), while insensitive to the spatial scales (L pc ) involved, and 
with a rather weak dependance on the pressure. However, for 
the dense gas-rich disks in ULIRGs the expected turbulent pres- 
sures can reach ~(l-3)x 10 7 Kcm -3 . Thus, for <r =1.2kms -1 
and turbulent-driving at the largest scales L pc =25-50 (disk scale- 
heights in local ULIRGs, DS98) it would be T kin ~(120-140) K. 
Atomic line cooling and gas-dust thermal interaction can lower 
these temperatures somewhat (e.g. Papadopoulos et al. 201 1) which 
nevertheless remain high enough to power the extreme-excitation 
global CO SLEDs is these systems. 

In reality both CR and turbulent heating will occur in the ISM 
of ULIRGs, powering high-excitation global CO SLEDs which 
would then be irreducible to ensembles of PDRs and/or SNR- 
shocked regions. This "irreducibility" has indeed been noted in 
the past for small molecular gas reservoirs (~(10 6 -10 7 ) M ) in 
the Galactic Center (Yusef-Zadeh et al. 2007) and the nucleus of 
NGC 253 (Bradford et al. 2003; Hailey-Dunsheath et al. 2008), yet 
in ULIRGs it will involve 3-4 orders of magnitude larger molecular 
gas masses. Unique observational signatures of such CR-dominated 
regions (CRDRs) or turbulenfly-heated regions (THRs) have been 
recently discussed in great detail (Papadopoulos 2010; Meijerink et 
al. 201 1; Bayet et al. 201 1) and most will be accessible to ALMA. 
Strong thermal decoupling of gas and dust with T kin 3>T dus t, a 
common aspect of both heating mechanisms, can be revealed via 
high-resolution mm/submm imaging of dust continuum and line 
imaging of dense gas tracers dense gas (e.g. CS, HCN). High-J CO 
and 13 CO line (e.g. J=6-5, 7-6) imaging will be particularly valu- 
able for distinguishing between the volumetrically heated CRDRs 
and THRs versus the surface-heated PDRs as those responsible 
for powering the observed CO SLEDs in (U)LIRGs (e.g. Hailey- 
Dunsheath et al. 2008). 



7.2 CO SLEDs and X co values as SF mode indicators in 
galaxies 

The low-efficiency (i.e. low Lji'/M(H2)) and spatially distributed 
star formation in isolated gas-rich disks, and the high-efficiency one 
occuring in the compact merger-driven starbursts in ULIRGs can be 
thought as two distinct SF modes. The latter accounts for only 1% 
of the total SFR density in the local Universe, but could be domi- 
nant in the distant Universe (e.g. Hughes et al. 1998). However the 
discovery of very gas-rich disks (f gas ~0.5-0.6) with extended low 
efficiency star-formation, and a space density 10-30 times higher 
than dusty mergers at high redshifts (Daddi et al. 2007), suggests 
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this mode as the prevelant one in the distant Universe. Recent stud- 
ies of such disks reported a low CO(3-2)/(2-l) ratio and a near 
Galactic X co factor, which were considered independent indicators 
of a disk-like low SFE mode (Dannerbauer et al. 2009; Daddi et 
al. 2010). Given the vast potential of molecular line observations in 
the upcoming era of ALMA, it is important to determine whether 
global CO SLEDs (by far the easiest to obtain for galaxies) and the 
corresponding X co values are indeed reliable SF mode indicators. 

The r 32 /2i distribution obtained for our sample (Figure 
11) makes obvious that the r 32 /2i recently measured in a dis- 
tant gas-rich disk (Dannerbauer et al. 2009) is within the range 
found for the vigorously star-forming LIRGs. Furthermore, as dis- 
cussed previously, uniquely attributing MW-type ISM conditions 
to subthermal low-J CO ratios is problematic without additional 
constraints (e.g. 13 CO lines). Finally, a near-Galactic X co fac- 
tor does not necessarily imply low-excitation SF-quiescent ISM. 
Indeed Galactic X co values can be associated with both low- 
excitation SF-quiescent (T b ,i-o~10K, n~500cm~ 3 ) as well as 
dense warm gas (Tb,i-o~60 K, n~10 4 cm~ 3 ) in star-forming sites 
since XcoOcy^n/Tb.i-o- This is actually one of the main reasons 
for the robustness of X co in estimating the mass of large molecu- 
lar cloud ensembles which inadvertendly encompass SF-quiescent 
as well as SF gas in spiral disks (e.g. Young & Scoville 1991). 
Finally, any AGN-powering of high-I CO lines in otherwise SF- 
quiescent galaxies, CR and/or turbulent heating, inject further un- 
certainties into any estimates of the dense/SF gas mass fraction 
f d =M(n>10 5 cm~ 3 )/Mtot(H 2 ) (another potent measure of the SF 
mode) using CO SLEDs. 

For a SF efficiency per dense molecular gas mass 
e» ig =L^ ) /M S F(H 2 ) that is a near-constant (~(250-500) L /M ), 
and set by the physics underlying SF feedback onto the ISM 
(e.g. Scoville 2004; Andrews & Thompson 2011), the f d be- 
comes an equivalent and much better measure of SF modes with 
/ d (mergers/starbursts)~(5- 10) x f d (disks/1 ow-SF-efficiency) (e.g. 
Solomon et al. 1992; Gao & Solomon 2004; Wu et al. 2005). The 
latter asumes the HCN/CO J=l ratio r H cN/co as a good proxy, 
which is also more practical to obtain (low frequency observations) 
than high- J CO SLEDs. An extensive body of data in the local Uni- 



verse allows both comparative studies of SF modes using "raw" 
Thcn/co ratios, and their "calibration" in terms of actual fd values 
(e.g. Papadopoulos et al. 2007; Gao & Solomon 2004). Sensitive 
CO and HCN J=l-0 and other types of molecular line observations 
that can discern SF modes of galaxies at high redshifts will be pos- 
sible in the future with MeerKAT, the SKA, and ALMA (Geach & 
Papadopoulos 2012). 



7.2. 1 BzK galaxies as SF disks: well-excited CO SLEDs beyond 
J=3-2? 



For the gas-rich disks of BzK galaxies at high redshifts: L^*' 



.(1- 



IR 

4)x 10^ L , which corresponds to M SF (H 2 )~(0.4-1.6)x 10 10 M 
of high-excitation gas (for e», g =25OM /L ). This is ~(5-13)% 
of the M to t(H 2 ) reported for these systems (Daddi et al. 2010) 
with ~10% being average. Local spirals on the other hand have 
M S F(H 2 )/M t ot(H 2 )~(3-10)% (e.g. IZw 1 in our sample), while 
in ULIRGs M S F(H 2 )/M tota i(H 2 )~(20^50)% (e.g. Solomon et 
al. 1992; Gao & Solomon 2004; Greve et al. 2009). For a two- 
phase decomposition of molecular line emission (a SF-quiescent 
and a SF-active phase) the global CO (T b (J+l-J))/(T b (l-0)} line 
ratio ((..} denote velocity/area averages) can be easily shown to be 



O) 

rj+i,J = r^ M 



1 + 



r j+i,j x; 

r .I + l,J X 



1 + 



1-fd 



(18) 



where (1) and (h) denote the quantities for the low and the high 
excitation phase. From COBE measurements (Fixsen et al. 1999) 
we obtain r^ ,=0.6, 0.28, 0.10, 0.046, <0.028, <0.017 for 
J+l=2,3,4,5,6,7 (upper limits are 3<r) for the the inner Galaxy 
which we use, along with Xco=5Xi, to characterize the low- 
excitation, SF-quiescent phase. For the high-excitation phase we 
use the CO (2-l)/(l-0) (—1.2-1.3) and CO/ 13 CO J=2-l (~10) ra- 
tios obtained for the SF "hot-spots" of Orion A and B (Sakamoto 
et al. 1994) as constraints on a LVG one-phase model. The dens- 
est and warmest phase compatible with them (T kin =(125-150) K, 
n~3xl0 5 cm~ 3 , and K vir =7) yields a CO SLED template with 
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Figure 11. The CO (3-2)/(2-l) brightness temperature distribution for our 
sample (shaded area: sources with CO source size ^14"), with the range 
corresponding to a measured value for a gas-rich disk in a BzK galaxy at 
z~1.5 (Dannerbauer et al. 2009). 



Figure 12. The CO SLED of a gas-rich disk with 10% of its molecular gas 
mass in a SF gas phase (typical for BzK galaxies) with the typical 30% 
dispersion expected from measurement errors (black line). The CO SLED 
of the inner Galaxy is marked by the red lines (see 7.2.1). 



ry^ ,=1.35, 1.33, 1.30, 1.27, 1.25, 1.22 for J+l=2,3,4,5,6,7, and a 

corresponding Xoo =2.2 Xi. Then Equation 18 yields rj+i,j=0.75, 
0.49, 0.34, 0.29, 0.27, 0.26 for J+l=2,3,4,5,6,7. The r 65 and r 76 
values assume the upper limits for and r^g, while setting the 
latter to zero yields r^O.25 and r^=0.24 as minimum values. 

Figure 12 shows the CO SLED for a fiducial gas-rich disk 
with 10% of its molecular gas mass in the SF phase along with the 
one for the inner Galaxy. For typical measurement uncertainties of 
cr(r)/r~0.25-0.30, only ratios rj+i,j with J+1^4 markedly devi- 
ate from a Galactic CO SLED while lower-J lines will remain com- 
patible with it. It is thus obvious that for the average SF gas mass 
fraction in BzK galaxies their CO SLEDs will be significantly more 
excited from that of the Galaxy beyond J =3-2. Future observations 
of higher-J CO lines in such systems are necessary to verify this. 



8 CONCLUSIONS, AND SOME IMPORTANT OPEN 
QUESTIONS 

We present multi-J 12 CO and 13 CO line data of total line lumi- 
nosities for a sample of IR-luminous galaxies (N=36) which in- 
cludes the most intense starbursts known in the local Universe. We 
then conducted a detailed literature search for all other IR-luminous 
galaxies with available total CO line luminosities up to at least 
J=3-2, the transition where the dense and warm star-forming gas 
phase starts dominating the CO SLEDs. The final sample contains 
70 galaxies with IR luminosities of L^~(10 10 -2x 10 12 ) Lq, mor- 
phologies ranging from isolated disks to strong mergers, and is an 
excellent resource for follow-up molecular line observations (e.g. 
HCN, CS, HCO + ) and mm/submm imaging of the CO or other 
molecular SLEDs with ALMA (~2/3 of the sample is accessible 
from Llano Chajnantor). In the present study it was used to ob- 
tain a better overal picture of the average physical conditions of the 
molecular gas in star-forming galaxies, our findings are as follows: 



(i) There is a surprisingly wide range of CO line excitation, from 
subthermally-excitated low-J transitions, to thermalized ones up to 
J=4— 3, 6-5. A positive K-correction is assured only up to the J=3- 
2, 4—3 lines, beyond which large variations of the emergent IR- 
normalized (Lj+i j /Lj^ ) line luminosities are expected. 

(ii) The highly excited global CO SLEDs are found solely in 
starburst galaxies with Lir.>10 Lq, with the exception of two 
found in galaxies with low or moderate star formation but host- 
ing powerful AGN. Such SLEDs imply tremendeous amounts of 
molecular gas mass (~(few)xlO 9 M ) in a dense (>10 4 cm -3 ) 
and very warm (Tkm-^100 K) state. 

(iii) The high densities and temperatures, as well as the strongly 
unbound motions (K v i r S>l) often found for the molecular gas of 
ULIRGs underlie global CO SLEDs which can remain well-excited 
and become partly optically thin up to high-J transitions. In the 
Galaxy such SLEDs are found only in strongly irradiated gas near 
HII regions and shocked gas in SNR-molecular cloud interfaces 
and involving only ~(l-5)% of the mass of individual Giant Molec- 
ular Clouds. In some ULIRGs of our sample the warm and dense 
gas phase can contain >50% of their total molecular gas mass. 

(iv) Highly supersonic turbulence and high CR energy densities, 
both permeating the ISM of ULIRGs, can volumetrically heat large 
amounts of dense gas to high temperatures, and easily account for 
their observed highly excited SLEDs. In the compact, and heav- 
ily dust-enshrouded ISM environments of ULIRGs such heating 
mechanisms may actually dominate over the classical mechanism 
via the far-UV/optical radiation fields and the photoelectric effect. 
The resulting highly excited global CO SLEDs will be irreducible 
to ensemble averages of Photon-dominated Regions (PDRs), an ex- 
citing possibility that will be further explored with Herschel obser- 
vations of the complete J-ladder above J=4-3 for several ULIRGs 
in our sample (Key project HerCULES). 

(v) As expected for a sample of IR-selected (and thus vigor- 
ously star-forming) galaxies, only few low-excitation CO SLEDs 
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Figure 10. The range of the normalized Lj^j j=Lj+i,j /Ljr CO SLEDs of the sample, along with a few key systems (lines beyond J=4-3 are interpolated 
using one-phase LVG models). 



are found. For only two of them a cold and gravitationally bound 
phase with moderate/low densities (typical of SF-quiescent molec- 
ular gas) emerges as the most likely ISM state, while the rest of the 
"cold" CO SLEDs may still belong to starburst systems. Only 13 CO 
and/or CO J+l >-J, J+l>3 lines can reduce this well-known degen- 
eracy, and as a result the low-J part of global CO SLEDs (J=l-0, 
2-1, 3-2) cannot be a unique indicator of the SF mode (i.e. merger- 
driven compact starburst versus isolated-disk low-efficiency star 
formation) in galaxies. 

(vi) The gas-rich disks recently discovered at high redshifts are 
examples of galaxies whose "cold" low-J CO SLEDs have been 
used to indicate Galactic-type low-excitation ISM, yet higher-J 
(J=4— 3 and higher) CO transitions are expected to be highly-excited 
and indicative of much more vigorous star formation environments 
than those found in local spirals. 

8.1 Important open questions 

If strong turbulence and/or high CR energy densities are respon- 
sible for the large amounts of very warm and dense gas found in 
some ULIRGs, they will also change the initial conditions of star- 
formation in such galaxies, set deep inside UV-shielded dense gas 
regions. The new conditions can lead towards a top-heavy stellar 
Initial Mass Function (IMF) (Hocuk & Spaans 2010; Papadopou- 
los et al. 201 1) with ground-breaking implications for the interpre- 



tation of IR/optical SEDs and HII region optical-IR lines of such 
systems in terms of star formation rates. Nevertheless the current 
data, while strongly suggestive of different dominant heating mech- 
anisms operating deep inside dense molecular clouds in ULIRGs, 
they cannot determine their state in detail. Recent work (Bayet et 
al. 2011; Meijerink et al. 2011) demonstrates this to be possible 
using the high-frequency coverage provided by Herschel, and the 
increased sensitivity in the mm/submm regime soon to be provided 
by ALMA. 



Finally, the wide range of average ISM conditions found in 
star-forming IR-luminous galaxies will strongly impact the so- 
called X co factor. Indeed it becomes rather hard to argue in fa- 
vor of one convenient, ULIRG-appropriate, X co factor given the 
diversity of CO SLEDs found even within the ULIRG class. The 
strong supersonic turbulence alone will "re-settle" large portions of 
molecular gas mass in such systems towards high density phases. 
These would be traceable only via high-J CO and heavy rotor (e.g. 
HCN, CS) transitions while low-J CO lines will be dominated by 
a warm, diffuse and unbound gas phase that contains little mass 
and has misleadingly "cold" low CO ratios. These issues are now 
addressed in detail in Paper II (Papadopoulos et al. 2012). 
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